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CHAPTER  I 


INTRODUCTION  AND  OVERVIEW 


1.0  Background 

Assessing  the  reliability  of  series  and  parallel  systems  in  terms  of  their  component 
reliabilities  is  an  important  practical  problem.  The  component  reliabilities  are  usually  assessed 
via  life  tests  conducted  under  controlled  test-bench  environments  which  are  generally  static 
over  time.  The  assessed  component  reliabilities  are  then  used  to  estimate  system  reliability, 
assuming  the  independence  of  component  lifelengths  and  ignoring  the  effects  of  the  operating 
environment  on  them.  Testing  the  entire  system  or  even  individual  components  under  actual 
operational  environments  is  rarely  feasible.  The  operating  environment  varies  from  one 
application  to  another  and  is  typically  harsher  or  gentler  than  the  test-bench  environment. 
Furthermore,  for  a  particular  application,  the  operating  environment  may  be  dynamic;  that  is, 
it  changes  over  time.  For  example,  the  failure  of  one  component  could  impose  an  additional 

1 


stress  on  the  surviving  components  and  it  facto  constitute  a  change  in  the  operating 
environment.  When  the  effect  of  the  operating  environment  on  the  system  is  ignored  and  the 
component  lifelengths  judged  independent,  the  calculated  assessments  of  system  reliability  fall 
short  of  actual  observation.  Such  experiences  have  cast  doubt  on  the  applicability  of  reliability 
theory  to  practical  problems. 

Multivariate  distributions  describing  dependent  lifelengths  of  biological  and  engineering 
systems  have  been  proposed  by  Freund  [1961],  Downton  [1970],  and  Marshall  and  Olkin  [1967]. 
More  recently,  Lindley  and  Singpurwalla  [1986]  -  henceforth  LS  -  considered  the  simple  case  of 
a  two  component  parallel  redundant  system  which  operates  in  a  static  environment  different 
from  the  test-bench  environment.  They  proposed  a  family  of  bivariate  distributions,  a  special 
case  of  which  is  the  bivariate  logistic  distribution  of  Malik  and  Abraham  [1973]. 
Notwithstanding  cumbersome  algebra  and  uninspiring  detail  (see  for  example  Nayak  [1987]  ), 
generalization  of  the  LS  model  to  the  multivariate  case  is  straightforward. 

Our  objective  here  is  to  expand  upon  the  basic  theme  of  LS  in  several  directions,  the  most 
important  being  a  consideration  of  dynamic  operating  environments.  To  facilitate  this  and  to 
establish  some  notation  and  terminology,  a  brief  review  of  the  basic  model  of  LS  is  necessary. 

Suppose  that  m  components  in  a  parallel  redundant  system  are  judged  to  have  exponential 
lifelengths  with  parameters  Alo,  A2#,  ...,  Amo,  where  the  A’s  pertain  to  test-bench  environments 
and  are  assumed  specified.  The  Alo’s  (or  more  generally  the  Alo(t)’s,  t>0,  if  the  lifelengths  arc 
judged  to  be  other  than  exponential)  will  be  referred  to  as  the  baseline  failure  rates  of  the 
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components.  Assume  that  the  effect  of  the  common  static  operating  environment  is  to 
modulate  (i.e.,  to  increase  or  decrease)  each  Al0,  i  =  0,1,. .  ,,m,  by  a  common  factor  r),  where  tj 
is  unknown;  suppose  that  our  uncertainty  about  rj  is  described  by  a  distribution  G.  The  notion 
that  the  operating  environment  -  a  covariate  in  the  biostatistical  literature  -  influences  the 
baseline  failure  rate  is  due  to  Cox  [1972];  ramifications  of  it  are  in  Kalbfleisch  and  Prentice 
[1980].  It  is  our  uncertainty  about  tj  that  induces  dependence  among  the  lifelengths  of  the 
components,  Tu  T2,  which  now  have  a  multivariate  distribution  whose  nature  is 

determined  by  the  assumed  form  for  G.  For  example,  when  G  is  taken  to  be  a 
distribution,  Tlt  T2l  have  the  multivariate  logistic  distribution  of  Malik  and  Abraham 

[1973]. 

The  behavior  of  the  reliability  of  series  and  parallel  systems,  under  the  set-up  of  LS  with  G 
a  gamma  distribution,  has  been  studied  by  Currit  and  Singpurwalla  [1987] /They  described 
circumstances  under  which  ignoring  the  influence  of  the  operating  environment  would  lead  to 
optimistic  or  pessimistic  assessments  of  system  reliability.  The  LS  model  gives  the  user  an 
example  of  the  importance  of  incorporating  dependencies  in  reliability  assessment  and  the 
consequences  of  a  naive  assumption  of  independence. 
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1.1  Extensions  of  the  Lindley-Singpurwalla  Theme 


Extentions  of  the  basic  theme  of  LS  can  occur  in  several  directions,  the  most  obvious  ones 

being  to  assume  that  the  baseline  failure  rates  are  no  longer  constant  and/or  to  choose  other 

forms  for  G,  say  the  lognormal  or  the  truncated  normal.  We  shall  not  pursue  these  lines  of 

development  because  the  underlying  analysis  is  straightforward  and  routine.  A  noteworthy 

extention  of  the  LS  theme  can  be  achieved  if  we  allow  the  operating  environment  to  be 

dynamic  by  making  17  a  function  of  time,  for  example  rj(t)  =  £  0k  t*  or 

*=o 

<  , 

i)(t)  —  exp  ^  1  }•  where  the  0ka  are  assumed  unknown.  In  the  former  case,  our 

*=o 

uncertainties  about  the  0’s  could  be  described  by  a  sequence  of  independent  gamma  or 
truncated  normal  distributions  whereas  in  the  latter  case  we  could  assume  normal 
distributions.  The  underlying  analysis  for  either  case  is  cumbersome  though  straightforward; 
consequently,  we  have  chosen  to  merely  state  the  ensuing  results  in  Appendix  B.  The  case 

I  . 

ij(t)  =  exp  {  53  Pk  *  }i  with  dependencies  between  the  0ks  (to  ensure  a  certain  type  of 
fc— 0 

smoothness),  has  recently  been  considered  by  Gamerman  and  West  [1987]  for  the  analysis  of 
medical  survival  data. 

Instead  of  assuming  that  q(t)  is  a  polynomial  in  t,  we  may  assume  that  rj(t)  is  a  piecewise 
constant  function  of  t  where  the  values  of  the  constants  are  unknown  and  are  described  by 
suitable  multivariate  distributions  which  ensure  some  form  of  dependence  between  the 
constants.  A  motivation  for  considering  the  piecewise  constant  function  is  its  meaningfulness  in 
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certain  applications  such  as  step-stress  testing  (see  for  example  Goel  and  DeGroot  [1979}  and 
Shaked  and  Singpurwalla  [1983]  ),  and  also  the  fact  that  certain  limiting  operations  on  it 
enables  us  to  describe  »j(t)  by  suitable  stochastic  processes;  as  we  shall  see,  these  processes  lead 
to  interesting  multivariate  distributions  for  the  lifelengths. 

Recent  work  has  explored  the  use  of  stochastic  processes  to  model  dynamic  environments. 
Flournoy  [1982]  has  considered  time-varying  covariate  processes  as  they  effect  both  failure  and 
censoring  mechanisms.  Cinlar  and  Ozekici  [1987]  have  used  a  semi-Markov  process  to  model 
the  environment  with  a  deterministic  function  relating  the  environment  to  the  failure  rate.  We 
will  focus  on  using  stochastic  processes  to  produce  useful  multivariate  distributions. 

In  what  follows,  we  shall  refer  to  ij(t),  t>0,  as  the  environmental  factor  function  (EFF). 
The  EFF  is  a  parameter  that  has  little  if  any  direct  physical  meaning.  It  is  introduced  for 

<r 

convenience  to  capture  our  opinions  about  the  physical  effects  of  the  operating  environment  on 
the  system.  However,  the  EFF  is  important,  since  the  nature  of  the  induced  component 
dependencies  depend  upon  the  form  of  the  EFF  and  the  manner  in  which  we  describe  our 
uncertainty  about  it. 

A  criticism  of  the  LS  model  is  that  the  operating  environment  is  assumed  to  have  the  same 
effect  upon  all  m  components  of  the  system.  This  assumption  is  manifested  by  our 
consideration  of  a  common  unknown  quantity  17  or  tj( t)  which  modulates  each  baseline  failure 
rate.  An  obvious  extension,  therefore,  is  to  assume  that  the  baseline  failure  rate  of  component  i 
is  modulated  by  an  unknown  quantity  (or  an  unknown  function  ij4(t)  ),  i  =  0,1,. ...m,  and 


5 


to  describe  our  uncertainty  about  the  q/s  at  any  time  t  by  a  meaningful  multivariate 
distribution.  This  idea  is  pursued  in  Chapter  2  for  the  static  environment  and  in  succeeding 
chapters  for  dynamic  environments. 


1.2  Outline  of  Dissertation  and  Overview  of  Main  Results 

In  Chapter  2  we  extend  the  LS  model  by  supposing  that  the  operating  environment  has 

different  effects  on  the  different  components.  To  do  this,  we  propose  a  multivariate  distribution 

which  describes  the  dependencies  between  the  EFFs  for  each  component  and  use  this  to 

produce  a  multivariate  distribution  for  the  lifelengths  of  the  various  components.  We  also 

✓ 

demonstrate  that  the  degree  of  dependence  between  the  EFFs  gets  reflected  as  inequalities  for 
the  survival  function  of  parallel  redundant  systems  (see  section  2.3).  The  entire  development  of 
Chapter  2  is  under  the  assumption  that  the  operating  environment  is  static. 

In  Chapters  3  and  4  we  consider  the  case  of  dynamic  environments,  with  Chapter  3  devoted 
to  the  case  of  the  EFF  described  by  a  piecewise  constant  function  and  Chapter  4  being  devoted 
to  the  case  of  the  EFF  described  by  a  continuous  time  stochastic  process  with  independent 
gamma  distributed  increments.  In  Chapter  3,  the  points  at  which  the  EFF  takes  jumps  may  be 
specified  or  random  and  the  sues  of  the  jumps  are  assumed  to  be  random  and  described  by 
gamma  distributions.  Considered  in  Chapter  3  (section  3.2)  is  the  case  in  which  the  EFF  is 
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described  by  a  “shot-noise”  process.  A  special  case  of  the  shot-noise  process  produces,  it 
appears,  a  new  family  of  multivariate  distributions  with  exponential  marginals.  In  section  3.3 
of  Chapter  3,  we  consider  the  case  of  a  dynamic  environment  described  by  piecewise  constant 
EFFs  which  has  different  effects  on  the  different  components  of  the  system. 

The  development  of  Chapter  4  produces  multivariate  distributions  which  appear  to  be  new. 
The  development  of  the  material  in  section  4.1  provides  an  alternative  derivation  of  the 
bivariate  exponential  distribution  of  Marshall  and  Olkin  [1967].  In  section  4.3  of  Chapter  4  we 
consider  the  case  of  a  dynamic  environment  described  by  a  gamma  stochastic  process  for  the 
EFF  which  has  different  effects  on  the  different  components  of  the  system. 
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CHAPTER  H 


DEPENDENT  LIFELENGTHS  INDUCED  BY  STATIC  ENVIRONMENTS 

HAVING  DIFFERENT  EFFECTS  ON  EACH  COMPONENT 

2.0  Introduction 

«■ 

The  model  proposed  by  LS  assumes  that  the  operating  environment  is  static  and  that  each 
baseline  failure  rate  Aiol  t  =  0,1, is  modulated  by  a  common  EFF  17,  with  the  uncertainty 
about  i)  described  by  a  specified  distribution  function  G.  A  natural  extention  of  the  above  is  to 
hypothesise  that  each  Aj#  is  modulated  by  a  different  EFF,  say  17  0  with  our  uncertainty  about 
17,-  described  by  a  specified  distribution  function  Git  t  =  0,1, It  is  desirable  to  assume  that 
the  m  environmental  factor  functions  are  dependent.  One  motivation  for  the  dependencies 
among  the  17/8  is  that  some  factors  which  constitute  the  environment  -  such  as  temperature  - 
may  have  an  identical  effect  on  all  components,  whereas  the  other  environmental  factors  have 
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different  effects  on  different  components.  A  plausible  model  for  describing  such  dependencies  is 
given  in  Section  2.1;  it  is  based  on  a  theme  proposed  by  Cherian  [1941]  and  David  and  Fix 
[1961]  -  see  also  Johnson  and  Kotz  [1972],  pp.  216-220. 


2.1  A  Model  for  Dependencies  Between  Component  EFFs 

To  simplify  matters,  we  shall  restrict  our  attention  to  the  case  of  a  two-component  system, 
so  that  >7i  and  tj7  are  the  only  EFFs  that  we  need  consider.  Assume  that  =  X„  +  X,,  where 
the  random  quantity  X0  captures  the  contribution  of  the  common  environmental  factors  on 
both  components,  and  the  random  quantity  X,-  captures  the  contribution  of  the  other  factors 
on  component  i,  i  ss  1,2.  The  Cherian- David- F ix  setup  assumes  that  X0,  Xx,  and  X3  are 
independent  gamma  distributed  variables  with  known  shape  parameters  a,  and  scale 
parameters  bi  for  i  =  0,1,2,  respectively.  Clearly,  rjt  and  rj7  are  dependent  with  a  joint  density 
at  Vit  >  0  given  by  the  Cherian- David-Fix  -  henceforth  CDF  -  bivariate  density. 

((,„,,)=  e'b| "bj ’>  ft  [  r-‘(«i)  ]  ft  b?‘ 

i  =  0  »  =  0 

min  (rji,fh) 

■f  (l[-*o)“rl  (li-Jo)”1*1  ,-(b»-b'-b>)'»  C,.  (2.1) 

xo=0 

In  what  follows,  and  throughout  this  dissertation,  we  shall  denote  the  gamma  distribution 
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by  writing  “X^  ~  G(  <*,,  b,-  ),"  where  the  gamma  pdf  is  of  the  form 


g  (  *  )  = 


_  b£ 


,a—\  —  b  x 


r(a) 


x  >  0  . 


2.2  The  Bivariate  Distribution  of  Lifelengtha 


In  this  section  we  again  focus  attention  on  the  two  component  case,  and  verify  via 


manipulations  identical  to  those  found  in  LS  that  under  the  Cherian-David-Fix  distribution  for 


the  EFFs  rjj  and  the  lifelengths  Tt  and  T3  of  a  two-component  parallel  redundant  system 
have  the  bivariate  survival  function 


F(  ri»  ri  )  =  p[  Ti  >  To  T2  >  r,  |  Xlo,  X7o ,  (  or^b,.  ),  <=0.1.2  )  = 


_ bp 

rl  +  ^2o 


r,  + 


rx  ,  rs  >  0  . 


(2.2) 


The  marginal  survival  function  of  T{,  <  =  1,2,  is 


?<(  Ti  )  =  Pt  T,-  >  rf-  |  Xiot  (  or,,  b,  ),  i=o,.  )  = 


_bo _ _~|ao  r  bi  >< 

~i  +  *>oJ  [A,.,  r,.  +  b J  ’ 


r4  >  0  .  (2.3) 
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Note  that  the  survival  function  (2.2)  is  the  product  of  three  contributing  terms  and  the 
marginals  (2.3)  are  each  the  product  of  two  terms,  one  representing  the  common  contribution 
of  X0  and  the  other  the  contribution  of  X{,  t  =  1,2.  The  development  of  (2.2)  and  (2.3)  has 
been  greatly  facilitated  by  the  assumption  that  X0,  Xt,  and  X3  are  mutually  independent.  If 
the  scale  factors  b0,  b1(  and  b3  are  equal  to  the  same  value,  say  b,  then  the  marginal  survival 
function  (2.3)  is  a  Pareto  with  parameters  (a0  +  or,-)  and  b/Alo,  i  =  1,2;  this  result  is 
analogous  to  that  of  LS. 

A  generalization  of  the  above  to  the  multivariate  case  is  straightforward,  although  it 
requires  the  specification  of  several  parameters. 

The  case  of  a  dynamic  environment  having  different  effects  on  each  component  is 
considered  in  Chapters  3  and  4. 


2.3  Inequalities  for  Survival  Functions  with  Increasing 
Degrees  of  Dependence 

An  examination  of  the  development  of  the  model  given  by  LS  and  the  model  (2.2)  given 
above  indicates  that  the  nature  of  dependence  between  Tj  and  T3  depends  upon  the  nature  of 
the  dependence  between  and  q3.  Note  that  in  the  LS  case  tfr  =  rj3  =  q,  so  that  the 
dependence  between  and  q3  is  the  strongest  possible,  being  one  to  one.  It  is  useful  to  view 
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dependence  between  two  variables  r/,  and  q2  subjectively,  in  the  sense  that  tj1  dependent  on  tj2 


implies  that  knowledge  about  causes  us  to  alter  our  assessment  of  uncertainty  about  q2 
made  prior  to  receiving  knowledge  about  rjj.  If  we  use  the  notation  X  =  Y  to  indicate  that  X 
has  the  same  distribution  as  Y,  a  simple  way  of  constructing  three  pairs  of  random  quantities 
(  Vit  »h  ),  (  *7j,  ),  and  (  r/t,  q2”  )  with  differing  degrees  of  dependence  is  to  introduce  four 

mutually  independent  random  quantities  X0,  Xlf  X0’,  and  Xj’,  with  X0(X1)  =  X^X^),  and 
define  =  q2  =  X0  +  Xu  jj2’  =  X0  +  Xj’,  and  ij2”  =  X0’  +  X^.  In  order  to  insure  that 
differences  between  survival  functions  are  caused  by  differing  degrees  of  dependence  and  not 
different  distributional  assumptions,  we  have  required  that  q2  =  tj2  =  q2”. 

If  we  assume  that  Xj  and  Xj’  ~  G(  ajt  bj  ),  «=o,i,  then  the  pair  of  EFFs  (  rju  q2  )  will 
result  in  the  bivariate  survival  function  given  by  LS,  where 


^Ls(  ri»ra  ) 


bj 

rl  +  *3o 


r„  r2  >  0  .  (2.5) 


If  we  assume  that  ax  =  ar2,  the  pair  (  q2’  )  will  lead  to  the  model  (2.2),  denoted  by 

f' CDp(  ri»  T a  )»  where  the  subscript  CDF  indicates  that  the  survival  function  is  based  on  the 
Cherian-  David-Fix  distribution  for  ijj  and  q2\  Since  and  i72”  are  independent,  the  pair 
(  lit  la”  )  wdl  lead  to  the  survival  function 
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find(  ’■i.'i  )  -  ;n  [bo  +\  r,]a° 


the  subscript  INO  denoting  the  fact  that  and  r;2”  are  mutually  independent.  Let 
(  Xx,  Yx  )  >  (  X2,  Y2  )  denote  the  fact  that  the  pair  (  X1(  Yx  )  is  more  (linearly)  dependent 
than  the  pair  (  X2,  Y2  ).  Then  (  «li»  *7a  )  >  (  »h*  »?»*)>(  ifu  »?a"  )•  It  is  easy  to  verify  that 
the  following  theorem  holds: 

Theorem  2.1 


FINd(  ti>  rJ  )  <  FCDf(  ri»  rJ  )  <  ^ls(  Ti>  r2  ) 
for  all  Tj  ,  r2  >  0,  and  Al0,  A2o  >  0. 

Thus  the  bivariate  survival  function  of  2-component  parallel  redundant  systems  increases  with 
the  degree  of  dependence  between  the  EFFs.  This  inequality  will  hold  for  the  multivariate 
survival  functions  of  all  m-component  parallel  redundant  systems  constructed  in  a  similar 
manner,  m  >  2.  If  we  set  =  r2  =  rm,  then  we  get  inequalities  for  the  system  survival 
function  of  m-component  series  systems. 
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CHAPTER  HI 


DEPENDENCE  INDUCED  BY  DYNAMIC  ENVIRONMENTS 

WITH  PIECEWISE  CONSTANT  ENVIRONMENTAL  FACTOR  FUNCTIONS 

3.0  Introduction 

In  this  chapter,  we  shall  assume  that  the  operating  environment  is  dynamic  and  that  the 
associated  EFF  is  described  by  a  piecewise  constant,  right  continuous  function  over  known  or 
unknown  time  intervals  [  t;,  t  .^.l  ),  j  =  0,1,...,  where  to=0.  Specifically,  we  let 
»?(*)  =  »  t  €  [t,-,  tj  +  l  ), 

with  the  unknown  and  assumed  to  have  various  distributional  forms.  Much  of  our 
development  will  be  based  on  the  assumption  that  the  dynamic  operating  environment  has  the 
same  effect  on  all  m  components.  We  shall  also  continue  to  assume  that  the  baseline  failure 
rates  of  all  m  components  are  constants,  so  that  for  t  6  [  t;,  t^,  ),  the  modulated  failure  rate 


of  the  »-th  component,  A,(t),  is  a  piecewise  constant  Xio  tjj.  Given  A{(t)  and ^  =  (  tlt  t2,  ...  ), 
the  lifelength  Tf  of  the  »-th  component  has  a  piecewise  exponential  distribution.  In  what 
follows,  we  shall  make  various  distributional  assumptions  on  the  ij^'s  and  _t ,  and  explore  the 
consequences  of  such  assumptions.  In  section  3.1,  we  assume  that  t,  is  known  whereas  in  section 
3.2  we  assume  that  it  is  not.  We  note  that  if  t.  is  known,  then  it  is  possible  to  determine  the 
intervals  in  which  any  two  times  and  r3  lie.  For  notational  convenience,  all  survival 

functions  are  stated  for  time  intervals  [  0,  rx  )  and  [  0,  r2  )  such  that 
0  <  tni  <  rt  <  t„i  +  |  <  t„,  <  r2  <  tnj  +  1. 


3.1  Describing  Each  as  a  Random  Sum  of  Innovations 

* 

As  a  starting  scenario,  we  assume  that  the  operating  environment  consists  of  at  most  5+1 

distinct  stresses  having  the  same  effect  on  all  m  components,  with  the  k-th  stress  contributing  a 

known  or  unknown  factor  (innovation)  ck  to  rjj.  We  also  assume  that  the  effects  of  the  factors 

ck,  k  =  0,1,... ,5,  are  additive  so  that 
5 

1j  =  Eb(*j)  c*»  j=o,i,... 

*= 0 

where  I*(t7)  =  1  if  the  fc-th  stress  is  present  during  ( ty,  ty_j_j  )  , 

0  otherwise. 

Let  pty  =  P(  I4(ty)  =  1  )  and  let  qky  =  1  —  pky.  Assume  that  p0j-  =  1,  V  j.  This 
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assumption  implies  that  factor  c0  is  always  present  in  the  operating  environment.  Note  that 
the  factors  ck  do  not  change  over  time;  it  is  presence  or  absence  of  a  factor  that  changes  the 
EFF. 

Were  we  to  assume  that  Ik(t and  ck  are  known  V  kj,  the  the  survival  function  of  Tt  and 
Tj  would  be 

**(  Tn  rj  )  =  P[Tt  >  rlt  Ta  >  r2  |  Alo,  AjCt  (  Ik(t^),  ck  ),  *=o,i . s,  j=o,i,...] 


-  n  /  ~  Iik^tn^  C*  +  +  -  Aj#  Ik(tn2)  C*  (Tj-tnj) 


e—  (A1«+A2«)  h(tj)  c*  (tj  +  i— li)  nrl  —  ck  (tj+i—tj)  1 

i=0  >=«!+!  / 


(3.1) 

s 

The  EFF  Vj  =  I*(t>)  ck  can  be  evaluated  for  each  time  period  I‘l'  ‘i  + 1>  when  the 

indicator  variables  Ik(t  ■)  and  the  factors  ck  are  assumed  unknown  by  first  conditioning  on  the 
ck's  known  and  then  averaging  over  their  distributions  Gk,  or  by  conditioning  on  the  It(ty)’s 
and  then  averaging  over  their  distributions. 

We  assume  in  this  analysis  that  the  innovations  are  mutually  independent;  that  is, 

iL  i»(‘,)  v  i,  **, 

i*(t i)  iL  v  *,  >*, 

c*  iL  c«  *9t» 
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where  “X  J]_  Y”  indicates  that  X  is  independent  of  Y. 


We  denote  the  vector  of  all  indicator  variables  as  J.  (t )  and  the  vector  of  all  factors  as  c. 


We  begin  by  conditioning  on  the  vector  c.  and  recognize  that  for  fixed  t,  j,  k, 
~  Ko  (  fcy)  Ck  j  _  J~~  ^  ti-l- 1  *f)  c* 


V,)1 * 


w.p.  ptj, 

w.p.  qtj. 


Since  our  indicator  variables  Ik(t^)  are  mutually  independent,  we  can  average  over  each  term 


to  arrive  at 


F(  Tj,t  j  |  £  )  =  P[Xj  Tj,  Tj  >  Tj  |  A]a,A]a,  c,  _t ,  p  ifj  j — o.i..  •  •*  * — o.i . a] 


■J30{[p‘ni 


~  [Alo(fl—  t»»i)  +  A2o(tni  +  1  trij)]  C* 


+  ‘itn,  1 


r  —  Ajo  (t2  —  tn3)  C4 

[  Pkn7  e  +  j 


ni  1  f  -  (Au+A,,)  (  t^  +  j-t;)  C* 

n  t  p*i «  +  j 

j=0 


"I1  .  -  Aj,  (  ti+l-»i)efc 


fl  t  P*ie 

J=nj+l 


+  1 


1 


(3.2) 


In  order  to  find  the  unconditional  bivariate  survival  function,  the  (n^-f-l)  terms  inside  the 
braces  must  be  evaluated.  This  would  be  cumbersome  for  large  itj  and  n2. 
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Alternatively,  we  can  condition  on  the  indicator  variables  and  write  the  conditional 


bivariate  survival  function  as: 

F(  rn  ra  I i(t)  )  =  P[ >  Tj,  T2  >  r2  |  Alo,  A2o,i(t), Ai  G*»  *=o,i,...  ]  = 


IfcCtnJ  tn1)  +  A2a(tfi^^.1— t„t)} 


Ifc(tf»2)  tn2) 


H-l 


""  1 


+  £  hftj)  (^lo  +  ^2o)(tj  +  i~tj)  +  12  !*(*>■)  +  *>)  1 

)— 0  j="l+l 


(3.3) 


where  Gk  is  the  Laplace-Stieltjes  transform  of  Gk.  The  bivariate  survival  function  is  not 
factorable  for  It(t;)  using  standard  nonnegative  distributions  Gt,  so  evaluating  the 
unconditional  bivariate  distribution  is  not  possible  except  under  special  circumstances  (such  as 
found  in  section  3.1.1).  Thus,  we  are  forced  to  consider  the  special  cases  presented  below.  In 
section  3.1.1,  we  examine  two  cases  that  lead  to  independent  rjj'a,  and  in  section  3.1.2,  we 
examine  two  cases  that  lead  to  dependence  between  the  ijS s. 


3.1.1  Independent  rjj'a  Described  as  a  Random  Sum  of  Innovations 

The  objective  of  the  first  part  of  this  section  is  to  produce  a  scenario  which  motivates  us  to 
describe  the  rj^'a  as  having  independent  gamma  distributions.  Such  a  setup  serves  as  a 
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foundation  foe  considering  the  gamma  process  (see  section  4.1.1)  to  describe  the  EFF. 

We  begin  by  assuming  as  before  that  the  indicator  variables  Ik(tj)  and  the  factors  ct  are 
unknown  and  are  mutually  independent.  The  distributions  will  be  stated  for  time  intervals 
[  0  ,  tv  )  and  [  0  ,  r2  )  such  that  0  <  tni  <  <  tni+1  <  t„,  <  r2  <  t  +l  ;  thus  we  are 

interested  in  the  EFFs  rj0,  ij2l  ....  rjn7.  We  assume  that  <  r,  this  prevents  us  from  using 
s+1  different  values  of  jj  to  solve  for  c0,  clf  ...  cs.  Therefore  P[  rjj  |  tj „,  ...,  ]  = 

P[  t)j  )  for  j  =  1,2,. ..,n2,  and  we  conclude  that  the  EFFs  are  mutually  independent. 

We  further  assume  that  the  number  of  stresses  present  during  time  interval  [  t;,  )  is 

known  -  denoted  by  Nj  -  although  the  identities  of  these  stresses  are  unknown.  If  we  describe 
our  uncertainty  about  each  factor  ck  with  an  identical  gamma  distribution  with  parameters  a 
and  b,  tj0,  t^,  . . rj„2  are  independent  gamma  distributed  random  variables  with  parameters 
N^a  and  b. 

Independence  between  the  rjy’s  may  not  reflect  the  actual  physical  system.  However,  this 
independent  model  serves  as  a  foundation  for  our  consideration  of  the  gamma  process  for  the 
EFF. 

The  development  of  the  multivariate  survival  function  follows  familiar  lines.  Specifically, 
for  the  case  m  =  2, 

rn  r7  )  =  P[  Ti  >  rj,  T2  >  r2  |  Al9,  A2o,  t.,  a,  b,  N;-  j=o,i . n2  ]  = 
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"  _ b _ 1N»ia  r _ b _ 1 

•*1  o(  Tl~ tf»l  )  +  A3o(  tni  +  l— tnl  )  +  b  (_  ASo(  TJ  —  t*»a  )  +  b  J 


lNnja 


"fiY. _ b _ iNJa  "ft1  r _ b _ iNJa 

i= 0  !_(Aio  +  A2o)(  S+i  “  )  +  b  J  j=n|+l  [_  Aa»(  fcj+i  ~  tj  )  +  b  J 


(3.4) 


The  marginal  survival  function  for  component  i,  t=l,2  is: 


F,-(  T,-  )  =  P(  T<  >  r4  |  Aio,  t,  a,  b,  Ni  j=o.i . ]  = 


r _ b _ -iN»<«  y  r  b  iN>a 

|_  A.«(  Ti~tni  )  +  b  J  j=0  Aio(  t;.+1  —  tj-  )  +  b  J 


0  <  tnj  <  r4  <  t„.  +  1 


(3.5) 


It  is  easy  to  verify  by  using  a  transformation  of  the  type  undertaken  by  LS  that  the 
marginal  distribution  function  of  either  component  reduces  to  a  piecewise  Pareto  distribution, 
so  that  (3.4)  may  be  referred  to  as  a  bivariate  piecewise  Pareto  survival  function.  The  survival 
functions  (3.4)  and  (3.5)  call  for  the  specification  of  too  many  hyperparameters  and  are  thus  of 
limited  practical  value  except  under  special  circumstances.  However,  as  stated  before,  in 
introducing  the  above  set-up,  we  lay  a  foundation  for  the  development  of  the  use  of  gamma 
processes  to  describe  the  EFF  -  such  processes  lead  to  interesting  results. 


Another  scenario  that  leads  to  independent  ijj’s  assumes  that  the  factors  ck  are  known  V  k, 
but  the  indicator  variables  It(t^)  are  not.  We  assume  as  before  that  the  I4(t^)’s  are  mutually 
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independent  for  k  =  0,1,... ,a,  j  =  0,1,...  Conditioning  on  c.  known  leads  to  the  survival 
function  (3.2).  Furthermore,  if  we  assume  that  ck  =  c,  k  >  0  and  pky  =  py,  k  >  0  (recall 
that  p0j  =  1  V  j  ),  then 

i  i*(t,)  ck  =  c  t  u(t>) 

k=0  k=0 

has  a  discrete  distribution  with  support  c,  2c,  ...,  (s+l)c  with  probabilities 

(  k  )  pi*  (1— pi)  S~k  >  k~  respectively. 


In  view  of  the  above. 


F(  ri»  ra  )  =  P[  Tj  >  rlt  T,  >  r2  |  Xu,  A2b,  t.,  c,  p,  j=o,i,...  ]  = 


t  ( • )  p„*  (i-p.,) — r  “+,)  [Ai-<r,'s>+w‘"+i“,"')l  c 


to{l)  P«a‘  (1-P»2)  a_‘ 


_fc  -  (‘+0  A2b  (ra— 

e 


r2 — t„2)  c 


a_k  -  (*+0  (^lo  +  ^Jo)  (  c 

)  e 


■fi‘  r  t  (j)p>  a--,) — 

i=»l+l  L  *=0  V  *  ' 


_k  ~  (‘  +  1)  A2o  (  tj.+l— 1;)  c 


(3.6) 


Alternatively,  if  we  assume  that  pt;-  =  pk  and  that  [  ty,  t^-  _j_  t  )  =  At  V  j,  then 
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F(  TU  T,  )  =  PfTj  >  Tlt  Tj  >  Ta  I  Xio,  X7o,c,  At,  pfc  k=0,i . a]  = 


n0  (i  p.  •" 


-tni)  +  A2o  At]  c*  +  qfc  J  (  Pk  e-  *U  (Tj  — tnj)  c*+<|t] 


[  P*  « 


~(^lo+^2»)  ^t  ct  ,  „  1  nl  r  ~  ^Jo  ^  ck  .  i  ni~ ("l  +  Ol 


+  q»]  1  [P»e 


+  <u>  ] 


Either  set  of  assumptions  reduces  the  number  of  parameters  that  need  to  be  specified  but  does 
not  significantly  reduce  the  algebraic  complexity. 


3.1.2  Dependent  q.-’s  Described  as  a  Random  Sum  of  Innovations 


The  two  special  cases  of  the  random  sum  of  innovations  scenario  contained  in  the  previous 
section  assumed  that  the  EFFs  tj0,  rj1,  ....  rjnj  are  mutually  independent.  A  more  realistic 
assumption  may  be  that  the  rjj'a  have  a  (time)  dependent  structure.  We  examine  two  cases 
below  that  lead  to  this. 

We  begin  the  first  case  by  assuming  that  the  indicator  variables  Ik(t^)  are  known  and  have 
the  following  special  structure: 


I*(t;)  =  0  0  <  >  <  * 


>=0.1,...  , 


1  j  >  * 


*  =  0,1...  .3  . 
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When  the  indicator  variables  have  this  structure,  the  EFF  rj-  for  the  time  period  [  ty,  t;.^_1  )  is 
3 

J  t  3  ~  ®ilf"  • 

4=0 

We  assume  that  each  factor  ck  is  unknown  and  that  ck  J]_  c,  k^i.  We  describe  our 
uncertainty  about  each  ck  through  a  specified  distribution  Gk.  Here  the  number  of  such 
distributions  is  equal  to  the  number  of  time  intervals.  A  set-up  like  this  would  be  relevant 
under  step-stess  testing  or  when  the  environment  is  described  by  a  stress  that  is  cumulative, 
such  as  doses  of  radiation.  Note  that  the  EFF  is  autocorrelated  in  time.  It  is  easy  to  verify  that 
cov[  q(s),  q(t)  ]  =  Var[  q(s)  ]  for  0  <  s  <  t. 

If  we  suppose  that  each  ck  is  distributed  as  G(  ak,  b4  ),  the  following  bivariate  survival 
function  is  obtained. 


?(  ri.  Ti  )  =  p[  Ti  >  rn  Tj  >  r,  |  Xl9,  Aj„I(t),i,  (  atj,  by  )  ;=o,i . n7  ]' 


”i  -  by  ”  ay  "a 

A  J*i  +  +  ^ao(r2—  *j)  ;  =  I|+1  J 


bj  +  ^2o(r2  —  t;) 


(3.8) 


This  special  case  can  be  transformed  into  a  shot  noise  model  if  we  assume  that  the 


transition  times  Tu  T},  . . .  can  be  described  by  a  Poisson  process.  For  details,  see  section  3.2. 


A  second  possible  set-up  assumes  that  each  »fy  is  dependent  upon  the  previous  *7y_t- 
Suppose  that  the  probability  that  any  given  stress  is  present  at  time  t  •  is  dependent  upon 
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whether  that  stress  was  present  at  time  ty  —  1.  The  stress  levels  remain  constant  (although 
unknown)  across  all  time  periods  so  dependence  will  arise  between  successive  time  periods  if  the 
indicator  variables  form  a  two-state  Markov  chain.  We  assume  that 


P[  =  •'  I  Ifc(to),  ]  =  P[  Ik(tj)  =  i  I  ]  j=  1,2,...,  IS  0,1;  and 

assume  the  following  transition  probabilities  for  each  /,  j  =  1,2,...  ,  and  for  fc=l,2,...,s  : 


P[I*(ti)  =  0|I*(ti_1)  =  0)  =  p(0,0) 

p[  4(‘i)  =  1 1  =  o  ]  =  p(o,i) 

P[Ik(ti)  =  0|Ik(ti_1)  =  l]=p(l,0) 
P[i*(ti)  =  i  \Wi-t)  =  1]  =  P(1,1) 


Note  that  p(0,0)  +  p(0,l)  =  1  and  p(l,0)  +  p(l,l)  =  1.  Thus  we  need  only  specify  two 
transition  probabilities.  The  one-step  transition  matrix  is  given  as  Pj,  below. 


PI  = 


for  i=t,a,... ;  *=o,i . s 


p(0,0)  P(0,1) 
p(i,o)  p(l,l) 

As  initial  conditions,  we  assume  that  the  probability  of  starting  in  state  0  is  equal  to  the 


probability  p(0,0)  and  that  the  probability  of  starting  in  state  1  is  equal  to  the  probability 
p(0,l)  for  k  =  1,2 ,...,».  We  continue  to  assume  that  the  factor  c0  is  always  present;  thus 

P[Io(t,)  =  1]  =  1  V;>0. 

To  solve  for  an  unconditional  survival  function,  we  note  that  for  any  time  period  j  and  any 
factor  c4, 

Ij(tj)  Cj  A,  (t.  i  ti) 

E[«  ’  ,+'  ’  =0|  = 

e  '*  A"  +  p<0,l)  +  1  ■  p(«,0)  ,  «n<i 
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—  4(^i)  ck  (^4.1 — ti) 

E[e  '  +  1  ,'|I*(ti_1)  =  l]  = 


-  ct  A,  (»i+1-ti) 


•  p(l,l)  +  1  •  p(i»o) 


Thus  for  any  j  >  1, 


E[e 

E[e 


U(kj)  c*  (tj-J-i—tj) 


|Ik(t,-i)  =  0] 

I  =  1  ] 


1 

0 

• 

1 

0 

e~  c*  A#  (tj+l-tj) 

1 

®ecau8e  [j  P(  !*(*(>)  )  —  ®  »  P(  I*(to)  )  =  l^]=[l,0^*Pj,we  can  evaluate  the  survival 
function  as  a  product  of  matrices  Pi;  as  defined  below.  Thus 


F(  rlt  r^sPlT^r,,  T,  >  r3  | 

Ai*»  ^a*»  ■£.*  -L>  p(0,l),  p(l,0)  ]  = 


—  c0  [  Alor!  +  A^r*] 


,ft{f.°i[|op.J][!]} 


(3.9) 


where 


q  e~ct(Alo+Aa#)(tJ.+|—  ty) 
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0 


;="i 

*=i . s 


1  o 

Q  e~c*(Aa*)(t,+i“ti) 


i  =  »j+l n2  — 1 

4  =  1 . 3 


1  0 

q  e~c»(^2o)(r2  —  ^nj) 


"  J 


4  =  1 . 5 


If  the  factors  cfc  k  =  0, are  also  unknown,  we  can  evaluate  the  unconditional  survival 


function  by  evaluating  the  product  of  the  matrices.  This  will  be  cumbersome  for  large  nlt  n2, 

at 

and  a. 


F(  ra  )  =  P[  Tj  >  rlf  Tj  >  r2 


^l*t  Gj»  * — 0,1 . s,  p(0,l),  p(l,0)  ]  — 


G'ol  Ai.rt  +  AJor,  ]  ^ft  |  Ec*  [  W  ]  [  .ft  P*i  ]  [  1  ]  } 


(3.10) 


where  Gq[s]  is  the  LST  of  the  distribution  of  c0. 


For  example,  suppose  0  =  t0  <  tj  <  rt  <  t3  <  r3  <  t3,  where  cA  ~  G(  ak,  b*  ): 
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■  h  {  P<M>*  +  P(0.D  P<1.0)  p(0,.)  [  b— (A,.+\.)(tl_to)] 

+  p(0,0)  P(0,1)  *1.0,  [^,.,^1;;+  >,.(1,_t|)]  •* 

+  rtu)  **  [  “ 

+  *0,0)  *0,1)  *1.1,  [  gjrw^^r,_tJ]  ”* 

+  p(0'1)  -X1'01 *«>  [  tnroa<^wr,-«j]  ** 

+  *1.0)  *!.!)■  [  B--rjr[— t  W+  A^-tj)  (3.H) 
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3.2  Modeling  the  EFF  a a  a  Shot- Noise  Process 


Suppose  that  the  operating  environment  consists  of  a  series  of  “shots”  (  or  “shocks”  )  of 
stress  whose  magnitudes  Xt,  k  =  0,1,...,  are  unknown.  Suppose  that  the  shocks  occur  over  time 
per  the  postulates  of  a  nonhomogeneous  Poisson  process  with  a  specified  rate  m(t),  t  >  0,  and 
that  whenever  a  shock  of  magnitude  X  occurs  at  epoch  s,  its  contribution  to  the  EFF  at  time 
s  +  t  is  X-h(t),  where  h(-)  is  a  positive  and  non-increasing  function  of  t  >  0;  h(-)  is  called  the 
attenuation  function.  Suppose  that  T^,  ...,  are  the  successive  epochs  at  which  shocks  of 

magnitude  XQ,  Xt,  . . .,  respectively,  occur;  then  the  EFF  at  time  t  is 

>7(0  =  2  X4  h(  t  —  T(k)  ) , 

where  h(u)  =  0,  for  all  u  <  0.  The  process  {  q(t),  t  >  0  }  is  called  a  shot-noise  process-,  see 
for  example,  Cox  and  Isham  [1980]  p.  135. 

If  the  baseline  failure  rate  of  each  component  is  constant,  then  modeling  the  EFF  as  a 
shot-noise  process  is  equivalent  to  modeling  the  failure  rate  of  each  component  as  a  shot-noise 
process. 

In  Figure  3.1  below,  we  show  the  sample  path  of  a  shot-noise  process  for  q(t);  to  ensure 
that  the  baseline  failure  rates  are  not  annihilated,  we  have  assumed  that  observation  on  the 
environmental  process  starts  at  T0  =  0,  at  which  time  the  first  shock  of  magnitude  X0 
necessarily  occurs. 
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EFF  ij(t) 


Figure  3.1  A  Sample  Path  of  a  Shot-Noise  Process 


The  use  of  a  shot-noise  process  for  modeling  life-lengths  is  made  possible  by  the  following 
two  assumptions: 

Al.  Xt  Jj_  T(0  ,  V  *  ;  this  assumption  states  that  the  magnitude  of 

each  shot  is  independent  of  when  the  shot  occurs. 

A2.  The  Xk’s  are  independent  and  identically  distributed  as  a  random 

variable,  say  X. 

Suppose  that  X  has  a  distribution  G,  and  that  G^(s)  is  the  Laplace-Stieltjes  transform  of 
G.  Then  the  following  theorem  specifies  the  survival  function,  in  the  operating  environment,  of 
a  single  component  system  with  a  constant  baseline  failure  rate  A„. 

Theorem  3.1 

Let  M(t)  =  f  m(u)  du  ,  H(t)  =  f  h(u)  du,  and  »;(t)  =  £  Xt  h(  t— T^j  ),  where 
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T(0)  =  0.  Then  under  Al,  A2,  the  lifelength  T  of  a  single  component  system  with  a  baseline 
failure  rate  X0  has  a  survival  function  for  r  >  0 


F(  r  |  A0  )  =  G£[  A,  H(r)  ]  e  M(r)  e 


|  G^[  A o  H(u)  ]  m(r— u)  du 


The  proof  of  theorem  3.1  is  facilitated  via  the  following  [  possibly  well  known  (  cf.  Lemoine 
and  Wenocur  [1986]  )  but  important  ]  lemma.  A  proof  is  provided  to  lay  the  necessary 
groundwork  for  the  proof  of  another  lemma  (Lemma  3.4). 


Lemma  3.2 

Suppose  that  T^j  <  T^j  <  •••  <  T^  are  the  epochs  of  n  occurrences  in  a 
nonhomogeneous  Poisson  process  with  a  rate  m(t),  t>0,  in  the  time  interval  [  0,  r  ).  Then  T,, 
i=  1, the  unordered  epochs,  are  independent  and  identically  distributed  as  a  random 
variable,  say  U,  whose  density  at  u,  0  <  u  <  r,  is  177-I,  where  M(r)  =  [  m(u)  du. 

M(r)  Jo 


Proof  of  Lemma  3.2: 

If  we  define  a  new  time  variable,  say  u,  where  U(t)  =  [  m(u)  du,  then  on  the  new  time  scale, 

J  0 

the  epochs  m(u)  du,  with  U^j  <  <  •••,  are  described  by  a  homogeneous 

Poisson  process  with  a  rate  1  (  see  Cox  and  Isham  [1980],  p.  48  ).  Therefore,  given  the  number 
of  shocks  n,  the  joint  distribution  of  the  ordered  epochs  <  U^2)  <  is  the  same 
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as  the  order  statistics  in  a  sample  of  size  n  from  a  uniform  distribution  on  [  0,  U(r)  ),  from 

which  it  follows  that  the  unordered  epochs  U{  ,  t  =  1,2, are  uniformly  and  independently 

distributed  on  [  0,  U(r)  ).  But  U(r)  =  M(r),  and  so  P[  U4  <  u  ]  =  P[  f  *  m(u)  du  <  u  ]  = 

Jo 

P[  M(T,)  <  u  ]  =  P[  T,-  <  M— 1(u)  ]  =  »  assuming  that  M— 1  exists.  It  now  follows 

that  since  P[  T,  <  M— *(u)  ]  =  —  ,  P[  Tj  <  u  ]  =  ,  and  thus  the  probability 

density  of  T4  at  any  u  is  .  □ 


Proof  of  Theorem  3.1: 

Let  A(r)  =  [  \0  ij(u)  du  and  suppose  that  we  condition  on  n  >  1  shocks  occurring  in  [  0,  r  ) 

Jo 

at  times  T  =  (  T(0)  <  T(l)  <  •  <  T(»)  )  with  each  shot  having  an  (unknown)  magnitude 

X;  recall  that  by  assumption  a  shot  must  necessarily  occur  at  =  0.  Then 

-  [a,  X  H(r)  +  A.  £  X  H(r-T(t))l 
F(r|A„,n,T,X)=e  L  *~l  J, 

which  because  of  Lemma  3.1  can  also  be  written  as 

fw.)  =  [«-a'xh(t)][«-  A-XH<r-u>]\ 

Note  that  the  above  step  would  not  have  been  possible  if  A2  did  not  hold.  If  we  remove  the 
conditioning  on  X,  then 

F(  r  |  A.,  n,  T  )  =  Gfc[  A0  H(r)  ]  [  G  A.  H(r-U)  ]  ]  "  , 
and  unconditioning  on  U,  we  have 

F(  r  |  A.,  n  )  =  G*x[  A«  H(r)  ]  £  £  G*x[  A„  H(r-u)  ]  jjjjj  du 
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Our  final  step  is  to  uneondition  on  n  >  0  -  which  has  probability  mass  function 
e  t.M(r.>  1  -  and  simplify  the  ensuing  expressions  to  get  the  desired  result.  □ 


The  function  h(u)  is  normally  chosen  to  be  a  decreasing  function  of  time  (  see  for  example 
Lemoine  and  Wenocur  [1985,1986]  )•  We  have  adapted  the  shot-noise  model  to  accomodate  a 
constant  stress  between  epochs,  consistent  with  the  postulates  of  this  chapter.  For  a  specific 
example,  suppose  that  X  ~  G(  a,  b  ),  h(u)  =  1,  and  m(u)  =  m,  for  all  u  >  0.  Then,  for  a 
£  1,  it  can  be  seen  that  for  r  >  0, 


» 


(3.12) 


and  for  a  =  1, 


mb 


-  1 


e 


mr 


Furthermore,  if  in  the  above  expression  we  set  m 


(3.13) 


F(  r  |  A.,  b  )  =  e~  mr  ,  (3.14) 

implying  that  the  lack  of  memory  property  of  the  component  lifelength  in  the  test  environment 
is  preserved  in  an  operating  environment  described  by  a  shot-noise  process,  if  the  shot 
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magnitudes  are  exponential  with  a  scale  parameter  b  and  the  shot  arrival  rate  is  the  ratio  of 
the  baseline  failure  rate  to  b.  This  result  is  intriguing  and  invites  some  form  of  interpretation; 
this  is  given  below. 


3.2.1  Behavior  of  the  Failure  Rate  Function 


We  can  easily  verify  that  the  failure  rate  of  (3.13)  is  for  r  >  0 


r(r)  = 


A o  (  1  +  mr  ) 
b  +  A«r 


^  (  1  +  mr  ) 


1 


4 .  *2  r 

+  b  T 


(3.15) 


thus  r(0)  = 


A* 

¥ 


and^lira^  r(r)  =  m.  Upon  examining  the  derivative  of  r(r), 


(  m  V’  —  I  ) 

namely  £  r(r)  =  g  ‘  T  ,  r  >  0, 

(y. +  T> 


(3.16) 


we  note  that  when  ^  <  (>)  m»  then  r(r)  >  (<)  0  for  r  €  [  0,oo  ),  thus  r(r)  increases 
(decreases)  from  ^  to  m;  when  ^  =  m,  ^  r(r)  =  0,  implying  that  r(r)  is  constant  at  m. 
Thus  the  behavior  of  the  failure  rate  function  depends  on  the  relationship  between  the 
parameters  ^  and  m;  see  Figure  3.2. 
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Failure  Rate  r(t) 


0  Time  r 

Figure  3.2.  Behavior  of  the  Single  Component  Failure  Rate  Function 

with  Varying  ^  . 

An  intuitive  explanation  of  the  monotone  increasing,  decreasing,  or  constant  behavior  of 
the  failure  rate  function  can  be  given  if  we  interpret  the  failure  rate  function  subjectively  -  as  it 
should  be  (  cf.  Singpurwalla  [1988]  ).  When  the  component  with  a  baseline  failure  rate  of  A0  is 
placed  into  the  operating  environment,  it  receives  a  shot  of  unknown  magnitude  X0;  this 
modulates  its  failure  rate  to  A«X0.  Since  X0  is  unknown,  it  is  reasonable  for  us  to  estimate  it  at 
its  expected  value  g,  making  r(0)  =  ^  .  As  time  passes,  we  learn  more  about  the  survival  (or 
failure)  of  the  item  and  this  plus  our  expectations  about  the  shot  arrival  process  causes  us  to 
revise  (upwards  or  downwards)  our  assessment  of  the  failure  rate  from  its  current  value.  For 
example,  if  the  operating  environment  is  such  that  there  are  no  shocks  on  (  0+,  oo  ) 
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[subsequent  to  the  first  shot  of  magnitude  X0  at  0  ],  then  the  expected  time  to  failure  of  the 
item  will  be  -p.  However,  if  shots  do  occur  (at  a  rate  m),  then  the  expected  time  between  the 
occurrence  of  the  shots  will  be  ^ .  At  the  occurrence  of  each  shot,  the  propensity  of  the  item  to 
fail  will  increase  and  so  will  our  assessment  of  the  failure  rate.  If  jp  =  i  equivalent  to  m  = 

we  will  expect  the  item  to  fail  due  to  the  effect  of  AoX0  at  about  the  same  time  as  the 
occurrence  of  the  first  shock.  Under  these  circumstances,  the  failure  rate  of  the  item  (as  viewed 
by  us)  will  be  constant  at  ^  because  (in  our  opinion)  the  dominant  failure  causing  mechanism 
will  be  due  to  the  A0  and  the  X0  combination  and  not  the  damage  due  to  the  shocks. 

If  if  <  m,  or  >  jg,  we  will  expect  the  first  shock  to  arrive  prior  to  failure  of  the  item 
due  to  the  effect  of  A0Xa;  given  that  the  item  has  not  failed  at  any  time  r,  our  assessment  of 
the  failure  rate  will  increase  since  we  expect  that  the  dominant  failure  mechanism  is  governed 
by  the  shot  noise  process.  It  is  helpful  to  conceptualize  the  situation  as  beinff  akin  to  that  of 
competing  risks  wherein  the  two  competing  mechanisms  of  failure  are  the  shot  process  over 
(  0^,  oo  ),  and  the  modulated  failure  process  determined  by  A»X0. 

If  ^  >  m,  or  ip  <  jg,  we  would  expect  the  first  shock  to  occur  after  the  item  has  failed 
due  to  the  effect  of  AoX0.  Given  that  the  item  survives  to  some  time  r,  our  judgement  is  that 
the  estimate  of  the  starting  value  of  the  failure  rate,  was  too  high  and  therefore  the  failure 
rate  should  be  lowered.  Note  that  (in  our  view)  the  dominant  failure  mechanism  here  is  not  the 
shot  noise  process.  The  arguement  used  here  is  analogous  to  that  used  by  Barlow  [1985]  to 


explain  why  mixtures  of  exponential  distributions  can  lead  to  distributions  with  decreasing 
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failure  rates. 


Finally,  under  both  cases  ^  >  m  and  ^  <  m,  the  question  arises  as  to  why  the  failure 
rate  will  asymptote  at  m.  We  can  explain  this  by  noting  that  in  a  shot  noise  process  with 
additive  damage  (like  the  one  considered  here),  the  time  to  failure  of  a  component  can  be 
viewed  as  the  time  at  which  the  cumulative  damage  exceeds  a  threshold.  Given  that  the  time 
between  the  shots  has  an  exponential  distribution  and  that  the  item  has  not  failed  at  some 
large  r,  its  failure  rate  at  r  will  behave  like  the  failure  rate  of  a  multiple  standby  system 
comprising  of  identical  components  with  exponential  lifelengths.  The  lifelength  of  such  a 
system  has  a  gamma  distribution  whose  failure  rate  asymptotes  at  the  failure  rate  of  each 
component  -  in  our  case,  the  quantity  m. 


3.2.2  The  Bivariate  Case 

The  bivariate  survival  function  of  a  two-component  parallel  redundent  system  with 
constant  baseline  failure  rates  Al4  and  A3o  is  provided  by  Theorem  3.3  below.  Its  proof  is 
facilitated  by  Lemma  3.4,  which  is  analogous  to  Lemma  3.2. 

Theorem  3.3 

If  M(t),  H(t),  and  Tj(t)  are  as  defined  in  Theorem  3.1,  then  under  A1  and  A2,  the  lifelengths  T, 
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and  Tj  of  a  two-component  parallel  redundent  system  have  the  bivariate  survival  function  for 


F(  I  *lot  ^2*  )  =  G£[  Au  HfrJ  +  A,.  H(r2)  ] 

[  1  G£  [  Xu  H(rj-nt)  +  A2o  H (r3-ut)  ]  m(ut)  d^ 

•  eJo 

fr2 

J  r  [  X3o  H(r2— u2)  ]  m(u2)  du2  _  M(Tj) 

•  e  *  •  e 


(3.17) 


Lemma  3.4 


Suppose  that  <  T^^jj  <  •••  <  T^njj  are  the  epochs  of  occurrences  in  a 

nonhomogeneous  Poisson  process  with  rate  m(t),  t  >  0,  in  the  time  interval  [  rlt  r2  ).  Then 


Tj  ,  «  =  (nj  +  1),  ...,  (n2),  the  unordered  epochs,  are  independent  and  identically  distributed  as 

"»(«) 


a  random  variable,  say  U,  whose  density  at  u,  rl  <  u  <  t2  is 


M(r2)— M(rx) 


Proof  of  Lemma  3.4: 

Following  the  proof  of  Lemma  3.2,  if  we  define  a  new  time  variable  U(t),  then  given  the 
number  of  epochs  in  [  rt,  r2  ),  the  joint  distribution  of  the  ordered  epochs 
^(nj  +  l)  ^  ^(nj+2)  <  <  ^(n2)  the  M,me  38  that  of  order  statistics  in  a  sample  of  size 

—  nt  from  a  uniform  distribution  on  (  U(rt),  U(r2)  ),  where  U(T,)  =  M(T,)  .  The 
remainder  of  the  arguement  is  the  same  as  that  of  Lemma  3.2.  □ 
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Proof  of  Theorem  3.3: 


Thin  proof  parallels  the  proof  of  Theorem  3.1  with  the  additional  proviso  that  for  Poisson 
proceaes,  the  number  of  events  in  [  rlt  r3  )  is  independent  of  the  number  of  events  in  [  0,  r,  )  - 
a  consequence  of  the  independent  increments  property.  □ 

Remark  3.5  The  marginal  distributions  of  T1  and  T3  follow  from  the  statement  of  Theorem 
3.3,  by  setting  either  or  r3  equal  to  0.  Specifically,  if  we  set  rx  =  0,  then,  for  r3  >  0, 


F(  r3  |  A3<>  )  —  G^[  X7e  H(t3)  ] 


* 


This  is  in  agreement  with  the  statement  of  Theorem  3.1. 

By  way  of  specifics,  suppose  that  X  ~  G(  1,  b  ),  h(u)  =  1,  and  m(u)  =  m,  for  all  u  >  0. 
Then  Theorem  3.3  leads  us  to  what  appears  to  be  a  new  family  of  bivariate  distributions  with 
exponential  marginals.  Specifically,  for  0  <  <  r3, 

F(  r|,fj  |  A1#,  A3#,  b,  m  )  = 


mb 


b  ~[  b  +  A3.  (r3-ri) 

_b  +  ^iori  +  ^2»r2_  _b  -1-  A1#rx  + 


b  +  A3„  (ra-rQ  ~|  (A„  +  A3o) 


[b  +  A3#(r3-r,)] 


.mb 

Jjj  —  m  r3 


(3.19) 


and 


(3.20) 


If  we  now  assume  that  Altf  =  A2o  =  A  and  that  m  =  g  ,  then,  for  0<rt  <  r2, 


?a(  ra  I  *a..  b,  m  )  =  [-+  b"~3] 


mb  _  i 

A28 


—  mr. 


Ti  > 


F(  tut7  |m)=, 


1  —  mr,  +  mr2 
1  +  mr,  +  mr. 


e 


m  r  2 


(3.21) 


and  for  ra  >  0, 

F2(r2)=  e“  mrj  ,  (3.22) 

an  exponential  distribution.  The  results  for  r2  <  r,  and  F  1(r1)  are  symmetric. 

The  results  of  this  section  can  be  generalized  to  multi-component  systems'  and  also  to  the 
case  in  which  the  occurrence  of  a  shot  at  T;  does  not  necessarily  lead  to  a  damage  Xj  ,  but 
associated  with  each  occurrence  of  a  shot  is  a  probability,  say  p,  of  a  damage  X j  ■  When  this  is 
the  case,  the  underlying  shot  occurrence  process  is  a  thinned  Poisson  process  and  the 
development  above  follows  mmtatia  mutandis. 
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3.3  Piecewise  Continuous  EFFa  Having  Different 
Effects  on  Each  Component 


In  this  section,  we  consider  an  explicit  dependent  structure  between  the  EFFs  for  each 
component  i,  *=l,2,...,m.  We  assume  that  each  EFF  ^(t)  is  piecewise  continuous,  sharing  the 
same  transition  times  .t  =  (  tt,  t3,  . . .  ).  The  Cherian-David-Fix  (CDF)  multivariate  model  is 
used  to  describe  the  dependence  between  components  during  each  period  [  t;,  tj  +  1  ).  For 
simplicity,  we  consider  only  the  two-component  parallel  structure.  Let  X0^  describe  the 
contribution  of  the  environmental  stress  common  to  both  components  during  the  time  interval 
[  t j,  tJ  +  1  ),  and  let  describe  the  contribution  of  the  other  stresses  to  component  t,  *=1,2, 
during  that  interval.  We  assume  that  the  X^’s  are  independently  distributed  as  G(  aijt  b^  ). 
Thus 


JL  x„ 

V  i,  and 

JL 

Vi,  i* ». 
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Gamma  Distributed  Variables 


If  we  assume  that  » =  X0j  +  XLj  and  tj7j  =  X0j  +  X2j,  the  bivariate  density  of 
(  li;*  *7a>  )  f°r  each  >  >  0  is  clearly  the  CDF  density  given  by  (2.1).  We  note  that  the 
component  EFFs  within  each  time  interval  are  dependent;  that  is,  rjl} r  and  rj7j  are  dependent 
for  each  j.  However,  the  EFFs  between  different  time  intervals  are  independent;  that  is, 
7tJ-  Jl_  *7/m  for  •,'=i,2  .  This  is  analogous  to  section  2.1  within  each  time  period 

( ti*  *>+ 1  )* 

The  survival  function  consists  of  the  product  of  contributions  of  the  form  of  (2.2)  for  each 
interval  [  t;-,  tj  +  1  ).  If  0  <  t1  <  r2  <  tx  ,  the  survival  function  is  identical  to  (2.2). 


F(  rnr2  )  =  P[  Tt  >  Tj,  Tj  >  t2  |  Alo,  Aj„,.t,  (  ai;-,  b^  )  <=0.1,2,  >=0,1,...  ] 


JOf»i 


•*lo(Tl  ~  ***1)  +  •^2o(tnl  + 1 —  *"l)  +  ^Onj 


bl"l 

Ol«l 

^2n2 

^l*(rl~  + 

^2o(t2  —  tnaJ  +  bjnj 
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("i-i  r 

{  n  ; 

t  J=0  [_< 


J0j 


(Al0+Aa#)  (t^j— tj)+b0j 


OtQj 


n 

1=0 


n2-l 

n 

;=0 


a7i  nrx 

•  n 

j*n!+l 


bQj  >Qi  \ 

Aj,  (tJ+i-t;)+boiJ  J 


(3.23) 


The  marginals  are  the  product  of  terms  similar  to  (2.3)  for  each  interval 
[  ty,  t i+l  ),  for  »  =  1,2: 


F t(rt)  =  P[  T,  >  r,  |  Afo,l,  (  or,,,  b„  )  /=o.,  i=o,i....  ]  = 


b0n,- 

a0ni 

b*ni 

A,o(r<— U.O  +  bon.^ 

Aj«(^i“t»j)  +  bjn^ 

v  /r _ hi_ 

i= 0  \  A i*  (4>  +  i 


)  +  b0i 


ao> 


_ IVj _  a»r  \ 

A««(tr+,  t>)+b,^  J 


(3.24) 


0  <  t„.  <  r,.  <  t„.+1 


If  the  scale  factors  for  all  component  distributions  are  equal,  that  is, 
b0;  =  bj;'  =  b3 j  =  bj  for  j  =  0,1,...  ,  then  the  marginal  distributions  (3.24)  are  piecewise 
Pareto  and  (3.23)  is  bivariate  piecewise  Pareto. 
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3.3.2  Component  EFFs  Modeled  as  Time-Dependent 
Gamma  Distributed  Variables 


If  we  assume  that  +  (  X0j-  +  X,;-  )  ,  i  =  1,2  ,  j  =  1,2,...,  with 

Vio  =  X0o  +  X,0,  where  each  Xi;  ~  G(  a^,  bi;-  )  and  all  X^-’s  are  mutually  independent, 
then 

j 

tfij  —  ^2  (  ^Ok  d"  ^ik  )  >  *  =  1.2  i  j  —  0.1 . 

k=0 

We  recognize  that  this  forms  a  time-dependent  model  similar  to  that  found  in  section  3.1.2, 
with  (  Tjjjt  tja  )  i  =  1,2,  autocorrelated  when  j^i  and  (  rjtj,  rj2j  )  dependent  for  each  j  >  0. 
The  bivariate  distribution  is  similar  to  (3.8),  with  separate  terms  for  X0j,  X1;  and  X2j- 


F(  rvr7  )  =  p[  Ti  >  Tv  T3  >  r2  |  Alc,  Aao,i,  (  aijt  b4>  )  .=0,1,2,  ;=o,i,...  ] 


n 

i—0 


' _ ^ _ >oi  f? 


bQj 


^3o  (r2~ tjj  +  boj 


*0  j 


nl 

n 

;=o 


bu 


^l#(Tl-*i)  +  blj 


n 

;=0 


^2*(r2~  t;)  +  bS> 


'2; 


(3.25) 
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3.4  Inequalities  for  Piecewise  Continuous  Survival  Functions 
with  Increasing  Degrees  of  Dependence 

We  can  examine  the  effect  of  increasing  degrees  of  dependence  between  components  1  and 
2  by  constructing  three  pairs  of  random  quantities  (  t)2j  ),  (  q1;-,  q2i’  ),  and  (  q^,  t ) 
for  each  j  =  0,1,  -.  in  the  same  manner  as  in  section  2.3.  We  define  mutually  independent 
random  variables  X0>,  X0>\  Xljt  and  Xlj*  with  Xoi  (Xiy)  i  X0/  (Xiy’)  and  define 
Hij  =  Hi]  =  X-oj  +  jt  Ojj*  =  4-  Xjj’,  and  q2j”  =  X0^’  4-  X j y  V  j  >  0.  We  note 

that  Tfaj  =  q 2j-*  =  q2j"  . 

If  we  assume  that  X,y,  XtJ’  ~  G(  a(j,  bfj  ),  i  =  0,1  ,  j  >  0,  then  the  sequence  of  EFFs 
{  (*?io»  »?2o  ),  (Oiu  «|ji  )i  ...  }  will  yield  the  LS  model  survival  function  (3.4),  with  a0j  + 

=  NjOt  and  b0j  =  bt>  =  b  V  j  >  0.  The  CDF  model  (3.25)  will  result  frorrf  the  sequence  of 
EFFs  {  (q10,  q20’  ),  (qa,  q2 ),...}  with  atl7  =  auj  and  b1;  =  b2j-  V  j  >  0.  The  sequence 
of  EFFs  {  (q10,  q20”  ),  (qa,  q21”  ),  ...  }  gives  an  independent  model  with  the  following 
survival  function: 


(3.26) 
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where  “IND”  indicates  that  71;-  Jj_  72j*  ,  j  >  0.  72y,  72j\  and  72y”  are  identically 
distributed  if  we  require  that  each  ari;-  and  b{j-  be  equal  across  all  three  models  (e.g., 
a00  (LS)  =  a00  (CDF)  =  ar00  (IND)  ).  Once  this  is  done,  it  is  clear  that 

(  7i jt  72 j  )  >  (  7i 72 /  )  >  (  7ij.  72 j”  )  J  —  Mv-,  and  that  the  following  theorem  holds: 

Theorem  3.6 

^INd(  rl»  Tl)  <  ^CDf(  rl*  r2  )  <  ^Ls(  rl>  r2  ) 
for  all  rt  ,  r2  >  0,  and  Al0,  A2o  >  0. 

Thus  the  bivariate  survival  function  of  2-component  parallel  redundant  systems  increases  with 
the  degree  of  dependence  between  dynamic  piecewise  continuous  EFF’s.  Theorem  2.1  is  a 
special  case  of  3.1,  when  rlt  r2  €  [  0,  tx  ).  This  again  generalizes  for  all  m-component  parallel 
redundant  systems,  m  >  2. 
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CHAPTER  IV 


DEPENDENCE  INDUCED  BY  DYNAMIC  ENVIRONMENTS  WITH 

ENVIRONMENTAL  FACTOR  FUNCTIONS  DESCRIBED 

BY  A  CONTINUOUS  TIME  STOCHASTIC  PROCESS 

4.0  Introduction 

In  this  chapter,  we  shall  continue  to  assume  that  the  operating  environment  is  dynamic 

and  that  our  uncertainty  about  the  EFF  ^(t)  or  its  cumulative  Y(t)  =  [  rj(u)  du,  is  described 

J  o 

by  a  continuous  time  stochastic  process  called  the  gamma  process.  The  gamma  process  for  the 
EFF  will  be  motivated  in  section  4.1.1  as  the  limit  of  a  piecewise  constant  EFF  (see  section 
3.1.1)  with  independent  gamma  distributed  innovations.  The  process,  being  nonnegative, 
nondecreasing  in  time,  and  possessing  independent  increments,  is  conceptually  appealing.  It  has 
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been  extensively  studied;  the  Gist  use  has  been  attributed  to  Moran  [1959].  In  this  dissertation, 
we  rely  primarily  on  results  of  Ferguson  and  Klass  [1972],  Cinlar  [1980]  and  Dykstra  &  Laud 
[1981];  the  use  of  gamma  processes  in  survival  analysis  is  primarily  due  to  Ferguson  [1973], 
Ferguson  Ic  Phadia  [1979],  and  Kalbfleiseh  [1978]. 

Definition  4.1  Let  a(t)  be  a  nondecreasing  left-continuous  real  valued  function  on  [0,oo) 
with  o(0)  =  0,  and  let  fi  6  (  0,  oo  ).  A  stochastic  process  {  Y(t),  t>0  }  is  said  to  be  a  gamma 
process  with  parameters  a(t)  and  /?,  denoted  “Y(t)  €  Gpr(  a(t),  P  )",  if: 

1.  Y(0)  =  0, 

2.  Y(t)  has  independent  increments,  and 

3.  Y(t)  —  Y(8)  ~  G(  a(t)  —  a(s)  ,  ^  )  for  any  0  <  s  <  t. 

The  distribution  G(  0,  )  is  defined  as  a  distribution  degenerate  at  0. 

Dykstra  and  Laud  [1981]  extend  the  gamma  process  to  include  a  time-varying  scale 
parameter  0(t). 

Definition  4.2  Let  /3(t),  t>0  be  a  positive  right-continuous  real  valued  function,  and  let 

Y(t)  €  Gpf(  a(t),  1  ).  The  process  Z(t)  =  i  /?(s)  dY(s)  is  an  extended  gamma  process 

1  o 

denoted  “Z(t)  €  Gpr(  o(t),  ^(t)  )." 
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Properties  of  the  gamma  procrag  (  Cinlar  [1980]  ;  Dykstra  tc  Laud  [1981]  ): 


Let  Y(t)  €  Gpr(  a(t),  0  ).  Then 

E[  Y(t)  ]  =  £  0  da(u)  and  Var[  Y(t)  ]  =  £  07  da(u) 


For  any  time  t,  Gy^»  the  Laplace-Stieltjes  transform  (LST)  of  the  distribution  of  Y(t),  is 


Gy(t)(s)  ~  [ 

GY(t)(8)  =  [ 


s  >  0  ;  with 


if  a(t)  is  continuously  differentiable. 


Properties  of  the  extended  gamma  process  (  Dykstra  tc  Laud  [1981]  ): 

Let  Z(t)  6  Gpr(  o(t),  0(t)  ).  Then 

E[  Z(t)  ]  =  £  0(u)  da(u)  ,  Var[  Z(t)  ]  =  £  0\u)  da(u),  and 

GZ(t)(8)  =  e*P  £  ~  £  ln(  1  +  80(u)  )  da(u)  J  >  8  >  0  , 

where  G^t)  “  ****  ^ST  of  the  distribution  of  Z(t). 


In  sections  4.1  and  4.3,  we  model  the  cumulative  environmental  factor  function 

Y(t)  =  f  rj(u)  du  as  a  gamma  process;  thus  the  cumulative  failure  rate 

Jo 

A(t)  =  |  A*(u)  dY(u)  is  a  gamma  process  if  Ae  is  constant,  and  an  extended  gamma  process 
if  A«  varies  with  time.  In  section  4.2,  we  model  the  EFF  itself  as  a  gamma  process. 
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4.1  Modeling  the  Cumulative  EFF  as  a  Gamma  Process 


We  begin  by  assuming  that  the  cumulative  EFF  at  time  r  is  Y(r)  €  Gpr(  or(r),  i  ),  with 

a(r)  continuously  differentiable,  and  ^a(r)  =  a (r).  The  cumulative  failure  rate  of  component 

i,  A i(r)  =  [  Aio(u)  dY(u),  t  =  1,2,  is  therefore  an  extended  gamma  process  with 

Jo 

parameters  (  a(r),  ). 

We  motivate  our  choice  of  the  gamma  process  for  Y(t)  in  section  4.1.1.  In  section  4.1.2,  we 
consider  the  special  case  of  the  baseline  failure  rate  Ajo(t)  being  a  constant.  This  leads  to 
several  bivariate  distributions  depending  upon  the  choice  of  shape  parameter  a(t);  a  special 
case  is  the  bivariate  exponential  distribution  of  Marshall  and  Olkin  [1967].  Thus  we  have  a  new 
derivation  of  this  well  known  distribution.  We  continue  by  considering  a  linearly  increasing 
baseline  failure  rate,  which  leads  to  a  (possibly  new)  family  of  multivariate  distributions  with 
marginals  that  have  a  logarithmic  failure  rate  function. 


4.1.1  Motivation  for  Modeling  the  Cumulative  EFF  a a  a  Gamma  Process 

In  section  3.1.1,  we  assumed  the  environmental  stress  to  be  constant  over  small  intervals 
[0,tx)  ,  [t1,t3)  ,  ...  .  We  captured  our  knowledge  about  the  stress  over  each  interval  [t;  ,  tj  +  1  ), 
j>0,  via  rjj.  Our  uncertainty  about  each  rjj  was  described  by  a  gamma  distribution  with 
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parameters  (  N;ar,  b  ),  with  tjj  independent  of  rj,  for  ail  j ^ «.  Thus  the  cumulative  failure  rate 
at  any  time  r,  0  <  t*»  <  r  <  tn_|_lt  for  the  sth  component  is: 

A.-(r)  =  \ie  £  rtj  ( ti+1  -  t,-  )  +  «;„(  r  -  t„  )  J  ,  i=i,a  . 

For  any  finite  positive  time  interval  At^  =  (  t  .^.1  —  t;  ),  it  is  possible  to  normalize  the  shape 
and  scale  parameters  of  the  gamma  distributions  at  each  time  t ^  by  defining 

Since  we  will  be  taking  the  limit  as  At  ^  tends  to  zero,  we  can  assume  that  the  intervals  are 
uniform  and  that  there  exists  an  integer  m  such  that  r  =  (m+1)  At.  Thus  the  cumulative 


failure  rate  can  be  written  as 
Aj(r)  =  A,- 


£  n j  At  ,  <= 

j— 0 


1,2. 


Since  rj ~  G(  a*(tJ)At  ,  b*At  ),  Aio  17^  At  ~  G(  a*(tj)At  ,  j-  )  and 


Aj(-r)  ~  G(  53  «*(tj)At  ,  y~  )  ,  with  characteristic  function 
j=o 


I — 7  .  -  £  ••(tjJAt 

^  _  /  ,  N-l  Aio  x  j=0 

#A,(r)  -  (  1  “  “b*“  } 


AsAtJ.0(mfoo),  the  characteristic  function  of  ^.(r)  ^>ecome8: 


4=1  X.  4T5a.  -j!  “» du 

-“a-  dAi(r>  =  o-V1  > 

Thus  for  any  r  >  0,  A4(r)  ~  G(  f  a*(u)  du,  £-  ).  If  we  consider  any  partition  of 

Jo  *i0 


[  0,  r  )  such  as  [  0,  s  ),  [  s,  r  ),  then  it  can  easily  be  shown  that 
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A  At)  —  \As)  ~  G(  f  a*(u)  du,  ),  with  A^r)  —  A,(s)  independent  of  A^s),  and 

J  8  Aio 

A{(0)  ~  G(  0,  y~  ),  implying  that  A<(0)  =  0.  Thus  we  may  conclude  that  the  cumulative 

•*to 

failure  rate  A,(r)  of  the  piecewise  constant  independent  innovations  model  of  section  3.1.1 
leads  to  a  gamma  process  as  the  uniform  time  intervals  (  t^ x  —  t;  ),  j  >  0,  tend  to  zero. 


4.1.2  Modeling  the  Cumulative  EFF  as  a  Gamma  Process  — 
The  Bivariate  Survival  Function  and  its  Marginals 


Recall  that  describing  Y(r)  by  a  gamma  process  with  parameters  a(r)  and  g  leads  to  an 

extended  gamma  process  for  A;(r)  with  parameters  <*(r)  and  *  =  l,2,...,m.  For  any 

* 

continuously  differentiable,  increasing,  and  real  valued  function  a(r)  and  any  continuous 
positive  real  valued  function  with  ^o(r)  =  a (r),  the  bivariate  survival  function  for 

0  <  r,  <  r2  will  be 


F(  Ti  )  =  P[  Ti  >  r„  Tj  >  r,  |  Ala(u),  A2.(u),  a(u),  b;  u  €  [  0,oo  )  ]  = 
exp{  -|T‘  ln(  1  +  Ale(uI+A3^u)  )  a(u)  du  } 

*  «P{  -|r|  *n(  1  +  )  a(u)  du  }  ;  (4.1) 

this  is  derived  directly  from  the  properties  of  the  extended  gamma  process. 

It  is  useful  to  examine  the  marginal  distributions  that  arise  from  this  formulation.  The 
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marginal  survival  functions  for  r {  >  0  arc 

p<(  Ti  )  =  pI  Tf  >  rj  |  Ajo(u),  a(u),  b;  u  €  [  0,oo  )  ]  = 

exp  {  —  |  '  ln(  1  +  — )  a(u)  du  }  ,  «  =  i,a  .  (4.2) 

Suppose  that  the  baseline  failure  rate  for  component  t,  Ajo(t)  =  Aio,  a  constant.  Then  the 
cumulative  failure  rate  A^r)  =  Ai#Y(r)  €  Gpr(  a(r),  ^  ),  i  =  1,2,  thus  for  0  <  rt  <  ra, 
p(  ru  t-j  )  =  p[  Ti  >  ri»  T,  >  r3  |  Alo,  X2e,  a(u),  b;  u  €  [  0,  oo  )  ]  = 


"  b  "I  «(ri)  f  b  I  “(ra)  “  «(ri) 
_  b  +  (Aio+Aj0)  _  L  b  +  . 


(4.4) 


with  marginal  survival  and  failure  rate  functions  for  r(-  >  0 
pi(  ^  )  S  p[  T,  >  t{  |  Aio,  a(u),  b;  u  €  [  0,oo  )  ]  = 


(4.5) 


ri(r.)  =  » (ri)  In  [  b  J  ,  *'=i.a  ,  respectively.  (4.6) 

The  marginal  failure  rates  consist  of  a  constant  times  the  derivative  of  the  shape  parameter 
a(r).  By  selecting  different  functional  forms  of  a(r),  each  suggested  by  different  physical 
models  for  the  environment,  we  can  derive  different  marginal  distributions  for  the  components. 
For  example,  if  we  describe  Y(t)  using  a  gamma  process,  we  might  motivate  it  by  describing 
the  EFFs  of  the  the  piecewise  constant  innovations  model  of  section  3.1.1  using  iid  gamma 
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variates.  When  we  derive  the  limiting  gamma  process,  the  shape  parameter  a(t)  will  be  linear. 

If  we  buppose  that  the  shape  parameter  is  linear,  say  a(r)  =  al  r,  ax  >  0,  then  the 
marginals  are  exponential,  and  the  bivariate  survival  function  for  0  <  rt  <  r3  is 
F(  t lt  r2  )  =  P[  Tj  >  Tj,  Tj  >  r7  |  Alo,  A2#,  atj,  b  ] 


f  b  -+•  X7o  f  b 

+j^0j  l  J 

b  +  A2o 

”iaiTt  r  ^  +  *i» 

>ir2  r 

>  +  Ai»+Aj0J  L*>  +  Ala  +  AjaJ  [ 

let  A *  =  at] 

i  r  ** + + 

■"L  •>  +  *>.  . 

V- 

*** 

III 

S* 

III 

0 

b" 

f  (b+Au)  (b+A2#)  ~] 
L  b  (b+Xt.+X79)  J 

,  then 

(4.7) 


P(  r,  )  =  exp  -  {  A?  Tj  +  X7  r3  +  X*7  max[  rltr7  ]  |  ,  tx,t7  >  0, 


(4-8) 


which  is  the  bivariate  exponential  (BVE)  distribution  of  Marshal]  and  Olkin  [1967],  with 


marginals 

Fj(r<)  =  P[Ti>r,|  Ajo,  ax,  b  ]  = 

«P  ~{  (  )  U  }  =  exp  -{  ax  In  [  b  \X*'  ]  r4  }  , 

r(  >  0,  « — 1,2  .  (4.9) 
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It  is  interesting  to  note  that  under  the  appropriate  assumptions,  modeling  the  cumulative 


failure  rate  as  a  gamma  process  can  lead  to  a  standard  shock  model  such  as  the  BVE. 

We  can  justify  this  development  as  follows.  We  know  that  the  gamma  process  is  a  pure 
jump  additive  process  -  therefore  if  A(t)  is  modeled  as  a  gamma  process,  we  can  represent  it  as 

OO 

A(t)  =  52  Xj,  where  X  -  is  a  jump  occurring  at  some  time  T;  €  [  0,  t  ).  Although  there  are  a 
i=i 

countably  infinite  number  of  jumps  occurring  in  [  0,  t  ),  if  we  let  N(  t,  e  )  denote  the  number 
of  jumps  in  [  0,  t  )  of  magnitude  >  e  for  any  r  >  0,  then  it  can  be  shown  (  see,  for  example, 
Basawa  and  Prakasa  Rao  [1980],  pp.  108-107)  that  {  N(  t,  e  ),  t  >  0  }  is  a  Poisson  process 
with  an  intensity  function  m(r).  If  we  consider  increments  of  size  ax,  then  jumps  of  magnitude 
X  >  kax,  k  =  1,  2,  ...,  arrive  in  accordance  with  a  Poisson  process  with  intensity  m(kax), 
and  we  can  proceed  by  decomposing  the  total  arrival  process  to  arrival  processes  of  jumps  of 
size  X  €  [  kax,  (k+l)ax  ),  which  will  have  intensities  am(  kax  )  =  m(  (k+l)ax  )  — 
m(  kax  ),  k  =  1,  2,  ...  . 

The  size  of  each  jump  is  a  random  variable  X  distributed  as  G^(x).  A  jump  in  the 
cumulative  failure  rate  represents  a  discrete  shock  to  the  system,  which  may  cause  the  failure 
of  one  or  both  of  the  components.  We  expect  that  the  probability  that  the  shock  causes  the 
failure  of  a  component  to  be  a  function  of  the  size  of  the  shock.  If  we  choose  an  increment  ax 
sufficiently  small,  we  can  regard  the  probability  of  failure  to  be  approximately  the  same  for  all 
X  €  [  kax,  (k  +  l)ax  ).  Let  p01(kax)  [  p10(kax)  ]  be  the  probability  that  a  shock  of  size 
X  €  [  kax,  (k  +  l)ax  )  will  cause  the  first  (second)  component  to  fail  and  the  second  (first)  to 
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survive.  Let  p00(kAx)  [  pu(kAx)  ]  be  the  probability  that  a  shock  will  cause  both  components 
to  fail  (survive).  Thus  we  have  for  each  arrival  process  of  shocks  of  magnitude 
X  €  [  kAX,  (k 4*1) ax  ),  k  =  1,  2,  a  standard  nonfatal  shock  process  with  (only)  a  common 
shock  applied  to  both  components.  Applying  familiar  techniques  (see  for  example  Mann, 
Schafer,  and  Singpurwalla  [1974],  pp.  148-149  ),  we  arrive  at  the  following  bivariate  survival 
function  for  ru  r3  >  0: 


PlTi  >  rlt  T3  >  r3  | X  €  [  knx,  (k4-l)Ax  )  ]  = 
exp{  -nm(kAx)  p01(kAx)  r,  -Am(kAx)  p10(kAx)  r3  -Am(kAx)  p00(kAx)  max(T1,r3)  }  , 


(4.10) 

which  represents  the  probability  that  the  components  will  survive  shocks  of  magnitude 
X  €  [  Icax,  (k  4-1)  ax  ).  The  probability  that  the  system  will  survive  all  shocks  is  the  product 
of  the  incremental  probabilities: 


F(  r7  )  =  P[  Tt  >  r„  Tj  >  r3  |  m(e),  p01(e),  Pi0(«),  P00(<);  e  >  0  ]  = 

f  oo  oo 

cxp^  -  £  ^m(kiix)  Poi(k^x)  rx  -  £  iim(kAx)  p10(kiix)  r7 
1  k=0  k=0 

oo  \ 

-  £  nm(knx)  p00(kAx)  max(r1,r3)  l  .  (4.11) 

t=o  J 

As  we  take  the  limit  as  ax  i  0,  we  have  for  ru  r3  >  0 

?(  tu  r3  )  =  exp{  -  |  p01(x)  dm(x)  r,  -  J  Pi„(x)  dm(x)  r3 

-  |  PooW  dm(x)  max(ri,r3)  }  .  (4.12) 


Thus  the  arrival  rates  in  (4.8)  are  equivalent  to  the  integral  terms  in  (4.12). 


To  summarize,  we  see  that  a  gamma  process  leads  to  the  Marshall-Olkin  bivariate 
exponential  distribution  derived  from  a  non-fatal  shock  model  with  only  a  common  shock 
applied  to  both  components.  The  major  difference  between  the  derivation  provided  by  Marshall 
and  Olkin  [1967]  and  this  derivation  lies  in  the  fact  that  the  original  derivation  assumed  that 
all  shocks  have  an  identical  probability  of  causing  failure,  while  we  assume  that  the  shock  (and 
thus  the  shock  failure  probability)  is  randomly  distributed. 


If  we  suppose  that  the  gamma  process  shape  parameter  at(r)  has  the  form  a(r)  =  ap  tp\ 
p,  ctp  >  0,  then  the  marginal  distributions  are  Weibull.  The  bivariate  distribution  for 
0  <  <  ra  corresponding  to  marginal  Weibulls  is: 

F(  T|,  ra  )  3  P[  Tj  >  Tj,  Ta  >  Tj  |  Aa#,  p,  ap,  b  ] 


b  +  Aa„  Hofprf 

b  larT a 

_  b  +  Alo  +  Aa„J 

_  b  +  Aao  J 

If  we  define  Af  =  ap  In  —  »  Aj  =  ar  In  ^  , 


and  A 


*  _  1  (b+Au)  (b+Aa>) 

ia  '  L  b(b+Au+A3.) . 


then 


?(  rlt  ra  )  =  exp  -  {  A?  rf  +  A}  r]  +  A?a  (max[  r„ra  ])'  }  , 

TvTi  >  0  »  (4.14) 
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which  is  a  bivariate  Weibull  distribution,  first  identified  by  Marshall  and  Olkin  [1967].  The 

marginal  survival  functions  are 

F|(  Ti  )  =  PI  T,  >  Ti  |  Aio,  P ,  af,  b  ]  = 


r,-  >0 


(4.15) 


Suppose  that  the  components  have  a  linearly  increasing  baseline  failure  rate,  denoted  by 

the  differentiable  function  Aio(r)  =  b  0iO  +  b  0tlt  for  some  b,  0iO,  0n  >  0,  i  =  1,2,  and 

suppose  that  the  cumulative  EFF  continues  to  be  described  1  y  a  gamma  process  with  a  linear, 

differentiable  shape  parameter  function  o(r)  =  a^r,  r  >  0.  Then  the  cumulative  failure  rate 

A ,(r)  =  [  Aio(u)  dYu  forms  an  extended  gamma  process  with  shape  function  a(r)  and  scale 
Jo 

A  r 

function  — — .  We  use  the  LST  of  the  increments  of  the  extended  gamma  process  to  compute 

the  survival  distributions.  For  0  <  rx  <  r2: 

rii  Ta  )  —  P[  Tj  >  rlt  Tj  >  r2  |  /?10,  0n,  07O,  /?21,  arj,  b  ]  = 


r  l±Ih+ hi  1  (I+ 

Ll  +  010  +  ^20  +  (011  +  02l)  rlJ 


n+02O  +  021rl~| 

^(i+a»)  r  i+#„+(J„r,  -I 

[_l  +  0JO  +  021r2 J 

l_l  +  01O  +  02O  +  (011  +  02l)  rlJ 

r  e 

-1  “l^a 

+  020  +  021 

(4.16) 
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This  appears  to  be  a  new  class  of  distributions  (at  least  in  a  reliability  context),  with  marginal 


survival  functions  for  >  0 


Fj(  Ti  )  =  p[  T<  >  Ti  I  0<O»  0il>  «!.•»]  = 


r  i±&.  1  ft ( 1+^'° >  r  „  i 

l_l  +  0iO  +  0ilr.J  + ^>0  +  ^.1  r.  J 

and  marginal  failure  rate  functions  for  rt  >  0 


(4.17) 


ri(r<)  =  <*i  In(  1  +  0.o  +  0ii  ri  )  .  <  =  »•».  (4-!8) 

which  are  logarithmic.  The  component  marginal  distributions  are  IFR,  with  the  rate  of  increase 
of  the  failure  rate  steadily  decreasing  over  time. 


4.2  Modeling  the  EFF  as  a  Gamma  Process 

In  the  previous  section,  we  described  the  cumulative  EFF  Y(r)  using  a  gamma  process.  In 
this  section,  we  describe  the  time- varying  EFF  rj(r)  itself  by  a  gamma  process  with  parameters 
(  a(r),  b  ).  This  section  is  motivated  by  the  work  of  Dykstra  and  Laud  [1981],  in  which  they 
model  the  component  failure  rate  using  a  gamma  process.  The  calculations  involved  in 
computing  the  characteristic  functions  of  distributions  with  stochastic  failure  rates  have  been 
developed  by  Antelman  and  Savage  [1965]. 

Although  the  gamma  process  has  independent  increments,  it  is  easy  to  show  that  the  EFF 
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»j(r)  is  autocorrelated  in  time.  Thus  Cov{  rj(s),  r/(t)  ]  =  Var[  >7(8)  ]. 

The  following  theorem  speciGes  the  survival  function  of  a  single  component  system  when 
the  EFF  is  described  by  a  gamma  process. 

Theorem  4.1.  (  Dykstra  &  Laud  [1981]  ): 

If  the  prior  over  hazard  [failure]  rates  is  Gpr(  <*(•),  /?( • )  ),  then  the  marginal  survival 
function  of  an  observation  X  is  given  by 

F(t)  =  exp  |  -  j  ln(  1  +  0(u)  (t-u)  )  da(u)  j  . 

If  arbitrary  limits  0  <  a  <  b  are  used  in  the  proof  of  theorem  4.1,  we  get  a  result  useful 
for  deriving  bivariate  distributions. 

Corollary  4.2. 

Let  Z  6  Gpr(  or(t),  0(t)  ).  Then 

Ez[  exp{  —s  |  Z(u)  du  }  ]  =  exp  {  —  J  ln(  1  +  s0(u)  (b— u)  )  da(u)  }  . 

In  section  4.2.1,  we  provide  a  motivation  for  using  a  gamma  process  to  describe  the  EFF. 
In  section  4.2.2,  we  present  general  results  for  bivariate  and  marginal  distributions  in  cases 
where  the  baseline  failure  rate  is  allowed  to  vary  over  time.  We  then  examine  a  special  case 
where  the  gamma  process  shape  parameter  is  linear  and  the  baseline  failure  rate  is  constant, 
which  provides  the  same  marginal  logarithmic  failure  rate  functions  as  found  in  section  4.1.2. 
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4.2.1  Motivation  for  Modeling  the  EFF  as  a  Gamma  Process 


In  Section  3.1.2,  we  considered  »j(t ;)  as  the  sum  of  j  +  1  independent  gamma  distributed 
random  variables  c0,  c j  with  parameters  atj  and  b.  It  is  possible  to  show  that  the  sum  ij(t) 
becomes  a  gamma  process  for  any  r  >  0  in  the  limit  as  the  time  interval  [  t;-,  tj+1  )  tends  to 
zero.  We  assume  that  the  intervals  [  tJt  t  )  =  At  V  j  >  0  and  that  there  exists  an  integer 
m  such  that  r  =  (  m+1  )  At. 

We  define  for  any  finite  positive  At  the  parameters 
and  b*  =  b  . 

Thus  Cj  ~  G(  a*(t;-)At  ,  b*  )  and  tj(t)  ~  G(  £  a*(t^)At  ,  b*  )  , 

)= 0 

r—  -  £  **(tj)At 

with  characteristic  function  =  (  1  —  )  1 


As  At  i  0  (and  m  T  oo),  the  characteristic  function  of  <f> 


becomes: 


lim 

m— *oo 


-£  «*(t;)At 

7=0 


Thus  for  any  r  >  0,  tj(r)  ~  G(  /J  cr*(u)  du,  b*  )  .  If  we  consider  any  partition  of 
[  0,  r  )  such  as  [  0,  s  ),  [  s,  r  ),  it  can  easily  be  shown  that 

tj(t)  —  »?(8)  ~  G(  Jg  ar*(u)  du,  b*  ),  ^(r)  —  »j(s)  is  independent  of  »?(s),  and  r;(0)  =  0. 
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From  the  above  argument,  we  conclude  that  the  EFF  » j(r)  of  the  cumulative  piecewise 
gamma  distribution  modeled  in  section  3.1.2  becomes  a  gamma  process  in  the  limit  as  the  time 
intervals  (  t;.  + 1  —  t;-  )  tend  to  wo. 


4.2.2  Modeling  the  EFF  as  a  Gamma  Process  —  The 
Bivariate  Survival  Function  and  its  Marginals 

The  EFF  r/(t)  continues  to  be  described  by  a  gamma  process  with  parameters  <*(t)  and  g. 
If  the  baseline  failure  rate  is  a  continuous,  positive,  real  valued  function  of  time,  then  we  need 
the  following  theorem  to  derive  the  bivariate  and  marginal  survival  functions. 


Theorem  4.3. 


Let  rj( t)  €  Gpr(  a(t),  g  ),  let  Ac(t)  be  a  known  continuous  positive  real  valued  function,  and 
let  A(t)  =  |  A«(u)  ij(u)  du.  Then 

F(  t  |  A0(t),  a(t),  b;  t>0  )  =  exp{  -  ln[  1  +  ^  J  A«(s)  ds  ]  da(u)  }  . 


Proof  of  Theorem  4.3: 


Partition  (  0,t  ]  into  n  x  ra  parts.  Let  h  =  g.  If  rj(t)  €  Gpr(  a(t),  g  ),  then 

r;(u)  =  |  g  dX(8),  where  X(t)  €  Gpf(  ar(t),  1  ).  We  can  approximate  the  cumulative 
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failure  rate  A(t)  by  a  double  sum;  that  is, 


A(t)  =  £  A„(u)  [  £  £  dX<»)  J  du  =  l±  A.(jh)  kg(  X(j£)  -  X(!i#)  ]  . 

^  Yi  —  x(l§)  ~  X(~1H~)  •  T^en  x,  -  G(  a(^)  -  a(TT“)t  1  )  .  and 

■  n  _  im 

MO  -  tE  A,,(jh)  £  Y,-  ,  and  we  can  expand  the  right-hand  side  as 

Dj=i  .=i 

5  {  [Mh) + Atf(2h)  +  •  •  •  +  A«(nh)]  £  Y,  +  [A„(2h)  +  •  •  -  A.(nh)]  £  Y,-  +  •  •  ■ 

* — 1  i—m  -f*  1 

+  [A«(nh)]  "g  Y<  } . 

«=(n  —  l)m  -f-1 

J  m 

We  introduce  for  convenience  a  term  W;,  where  W j  =  JT  Yf  ,  j  =  l,...,n.  Then 

*=(/  — l)m  +  l 

W j  ~  G(  a(jh)  —  a((j-l)h)  ,  1  )  and  Wj  JJ_  W4,  from  the  independence  of 

increments  of  the  original  gamma  process.  We  can  now  write  A(t)  as 
A(t)  £.£  i  t  ^(kh)  h  J  Wj  . 

Using  this,  we  find  the  LST  of  A(t). 


E[  exp{  -  A(t)  }  ]  ~  ft  E[  exp{-  £  [  £  A,(kh)  h  J  W ,  }  ]  .  Since 

j=i  °  k=j 

W,  ~  G(  aQh)  -  a((j-l)h)  ,  1  ), 

1  =  0  +  0 


E[  *  wi  1  _  ,  ,  ^  C  «0h)  -  «.(0-!)b)  } 


.  Thus 


E[  e  A(t)  1  =  ft  1  1  +  i  £  A.(kh)  h  ]  1  ~  1  .  To  got  on  oxnct  volo.  fo, 

E[  e  ^  ],  we  take  the  limit  of  the  approximating  sums  as  a  |  oo  (  thus  h  [  0  ).  Therefore 

E[  «" A(t)  I  =  Km  exp{  -  £  [  aQh)  -  a(G-l)h)  ]  ln[  1  +  I  £  A„(kh)  h  ]  }, 

>= i  u  *=; 

and  the  survival  function  is 


j- 1 


F(  t  )  —  E[  e  ^  ]  =  exp{  —  J  In[  1  +  g  J*  A»(s)  ds  ]  da(u)  } 


0 
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Note  that  if  A0(t)  =  A«  V  t  >  0  (constant),  then  Theorem  4.3  is  equivalent  to  Theorem 
4.1. 

The  bivariate  survival  function  for  0  <  rt  <  r3  follows  directly  from  Theorem  4.3.  Thus 
F(  r„  r3  )  =  P[  Ti  >  ru  T,  >  r3  |  Ala(u),  A3o(u),  a(u),  b;  u  €  [  0,oo  )  ]  = 

F(  ,r3  )  =  exp{  -  }Q  1  Inf  1  +  b  {y1  (Aio(®)+A2»(8))  *  1  *<u)  du  ) 

•  exp{  -  |r|  N  1  +  g  Ju2  A3»  ds  ]  a(u)  du  }  ,  (4.19) 

with  marginal  survival  functions  for  rf  >  0 

F(  )  =  P[  T,-  >  r,-  |  Aia(u),  a(u),  b;  u  €  [  0,oo  )  ]  = 

«p{  ~  |o’  M  1  +  g  |u*  ds )  a(u)  du  }  , 

and  failure  rate  functions  for  r {  >  0  (4.20) 

*(*■,)  =  \Tq{  ( I  +  l  As»  ds  J~J  a(u)  du  ,  ,  =  !.:»  .  (4.21) 

Once  again,  we  can  chooee  plausible  functional  forms  for  a(t)  and  Aio(t),  suggested  by  the 
physical  model  of  the  environment,  that  allow  for  a  closed-form  solution  for  the  distribution 
functions.  Choices  of  a(t)  said  Ajo(t)  that  lead  to  interesting  marginal  distributions  are 
provided  below,  along  with  their  respective  bivariate  distributions. 

One  special  case  of  interest  in  our  motivation  of  the  gamma  process  (section  4.2.1)  is  the 
process  as  a  limit  of  the  sum  of  iid  gamma  variates.  If  c;  is  identically  distributed  for  all  j, 
j  =  0,1,...  ,  with  »j(t)  =  53  cj  >  then  rj(t)  —  rj( s)  =  53  cj  be  described  by  a  gamma 
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distribution  with  a  shape  parameter  that  depends  on  the  length  of  the  interval  (t-s),  which  will 
lead  in  the  limit  to  a  gamma  process  with  a  linear  shape  function  a(t),  say  a(t)  =  t,  for 
some  ot!  >  0,  t  >  0. 

We  assume  that  r/(r)  €  Gpr(  a,r,  g  ),  which  implies  that  the  component  failure  rate  is 
A^r)  =  Aio  ij(r)  €  Gpr(  alrl  ^  ).  For  convenience  in  notation,  we  let  ^  and 


The  bivariate  survival  function  for  0  <  rt  <  r7  can  be  found  using  Corollary  4.2. 


F(  T’l*  r2  )  =  P[  >  rlt  Ta  >  r2  |  f3n,  P71,  <*i  ]  — 


—  (  1  +  (0U  +  02l)Tl  ) 


(0n+0ai) 


_  ai 

(  1  +  /?«(  ra-rj  )  ) 


’  1  +  02l(  r2_ T1  )  1  <*lrI 

r  e  1  aiTi 

_  1  +  (0U  +  /*2l)  rl  _ 

L  1  +  ^ai(  Ta  — ri  )  J 

(4.22) 


The  marginal  survival  functions  for  rt  >  0  are 


F(  rf)  =  PlT^rJ^j,  ax  ] 


Qi 

SI 


=  < 1  +  [irfcrr] air< 


(4.23) 


with  failure  rate  functions  for  r,  >  0 
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r,(rJ  =  <*1  In[  1  +  PaTi  ]  ,  i= 1,2  . 


(4.24) 


It  is  interesting  to  note  that  the  marginal  distributions  are  the  same  as  those  obtained  in 

the  previous  section  using  an  extended  gamma  process  with  shape  parameter  a(r)  =  atr  and 
A-  (r) 

scale  parameter  =  fJix  r.  However,  we  have  a  different  bivariate  distribution  leading  to 

this  type  of  marginal  distribution. 

If  we  suppose  that  the  component  baseline  times- to-failure  are  described  by  independent 
Weibull  distributions,  then  Al0(t)  =  p  Af„  tf  1  p> o  *  =  1,3.  The  marginal  failure  rate 
functions  have  closed  form  expressions  if  the  shape  parameter  er(t)  of  the  EFF  gamma  process 
is  a  p— degree  function  of  time.  If  we  suppose  that  a(t)  =  ap  tp  ,  a,  >  0  ,  then  a(t)  =  ^o(t) 
=  p  atp  tp  \  and  the  failure  rate  is 

r<(t)  =  p  ap  tp  1  ln[  I  +  tp  ]  ,,  ap  >  0  ;  t  >  0,  <= 1,3  .  (4.25) 

For  example,  if  we  let  p  =  1,  then  o(t)  =  ajt,  a(t)  =  alf  and  Afa  t*  =  Aja  t.  Thus 
ri(t)  =  or1ln(I  +  ^t],  <  =  1.2, 

which  is  consistent  with  the  results  obtained  previously.  Bivariate  distributions  can  easily  be 
derived  using  Theorem  4.3. 
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4.3  Modeling  the  Cumulative  EFF  as  a  Gamma  Process 


Having  Different  Effects  on  Each  Component 

In  this  section,  we  use  the  approach  embodied  in  the  Cherian-David-Fix  class  of  bivariate 

distributions  to  define  a  dependent  set  of  gamma  processes,  assuming  that  some  environmental 

stress  elements  affect  the  components  in  the  same  manner,  and  other  stress  elements  affect  the 

components  differently.  We  describe  the  cumulative  effect  of  the  common  stress  at  time  r  with 

a  gamma  process  denoted  as  Y0(r);  that  is,  Y0(r)  €  Gpr(  a0(t),  g  ).  We  describe  the 

cumulative  effect  of  the  stress  unique  to  component  i  at  time  r  with  a  gamma  process  denoted 

as  Yj(r);  that  is,  Y{(r)  €  Gpr(  ait  g  ),  i  =  1,2.  We  also  assume  that  the  processes  have 

mutually  independent  increments;  thus  for  any  0  <  s  <  t  and  any  0  <  u  <  v, 

Yy(t)— Yy(s)  J|_  Yt(v)— Y,(u),  >^i  ,  >,>=0,1,2.  The  component  baseline  failure  rates  are 

assumed  constant,  so  A,(t)  =  [  Aio(u)  ^(u)  du  =  Ajo  [  Y0(t)  +  Yj(t)  ],  «'  =  1,2.  This  model 

Jo 

yields  the  bivariate  survival  function  for  0  <  fj  <  ra  as 

F(  ri>  Ti  )  =  p(  Ti  >  ri>  T,  >  r,  |  Aj,,  \Jo,  Qfi(u),  b;  u  €  [0,oo),  .=0,1,2  ] 


r  b 

«o(ri)  f  b  1*o(*’a)-«o(ri) 

r  b  i«i(ri) 

_b  +  Alo  +  A2o_ 

|_b  +  A2#_ 

b  +  Ala 

(4.26) 
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with  marginal  survival  functions  for  r{  >  0 

F(  Tf  )  =  P[  T,  >  rt  |  Xit ,  or,(u),  b;  u  €  [0,oo),  i=o,i  ]  = 

T  b  lao(ri)+a,(ri) 

-  L^ns-J 

and  marginal  failure  rate  functions 

r<(ri)  =  (  *o(ri)  +  »i(r<)  )  In  >  *«>.*. 

The  marginal  results  obtained  here  are  essentially  the  same  as  those  obtained  from  Section 
4.1.2.  The  marginal  distributions  may  be  IFR,  DFR,  or  have  mixed  increasing  /  decreasing 
failure  rates,  depending  upon  the  shape  functions  a^r),  t  =  0,1,2,  chosen.  For  example,  if  the 
processes  Yj(r)  have  a  linear  shape  function  a{(r),  the  marginal  distributions  are  exponential, 
and  if  the  shape  functions  are  of  the  form  aip  tp,  aip  >  0,  p  >  0,  then  the  marginal 
distributions  are  Weibull. 

4.4  Inequalities  for  Survival  Functions  of  Continuous  EFFs 
With  Increasing  Degrees  of  Dependence 

The  bivariate  survival  functions  can  be  compared  under  increasing  degrees  of  dependence 
for  the  continuous  distributions  of  sections  4.1  and  4.3.  Let  Y0(r),  Y0’(r),  Yj(r),  Yt’(r) 
denote  gamma  processes,  with  Y 0(r)  ^  Y0’(r),  Y^r)  ^  Yi’(r)i  and  Y^r),  Y/(r)  mutually 


(4.27) 


(4.28) 
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independent,  i  =  0,1,  V  r  >  0.  We  assume  Y,(t)  €  Gpr(  ar,(t),  g  )  and  we  construct  three 

pain  of  random  functions  (  A1(t),  A2(t)  ),  (  A^t),  A2’(t)  ),  and  (  At(t),  A2”(t)  )  by  defining 

^  =  Y0(t)  +  Y1(t),  =  Y0(t)  +  Yj’(t),  and  ^  =  Y0’(t)  +  Y1»(t). 

In  this  section,  it  is  more  convenient  to  work  with  the  cumulative  failure  rate  functions  A(t) 

ft  A  ft) 

than  the  EFFs  r/(t),  noting  that  the  cumulative  EFF  for  component  t  is  0i(u)  du  =  -rv  , 

J  0  Aio 

i  =  1,2. 

The  LS  model,  where  t  yields  the  survival  function  for  0  <  <  r2 

*10  A7o 

PLs(  ri»r3  )=  p[  Ti  >  T3  >  r3  I  Xio  X7o,  “o(“)>  ai(u).  b;  u  €  [0,oo)  ] 


_  r  b 

_b  +  Ala  +  A2o_ 


ao(Ti)+«i(ri) 


f  b  l[oo(fa)-bat(f2)]-[a0(r1)-ftti(ri)]  (4  2g) 

|_b  +  A2oJ 

The  Cherian-David-Fix  (CDF)  model  using  At(t),  A2’(t)  yields  for  0  <  <  r2 

PCDF(  ri>r3  )=  p[  Ti  >  ru  T2  >  r2  |  Al8,  A2„,  o0(u),  <*i(u),  b;  u  6  [0,oo)  ]  = 


r__b__l0ro(n)  r  b_l(“o(r3)+ai(r2)]-Oo(ri)  [  blai(Ti) 

[_b  +  Al8  +  A2oJ  j_b  +  A3oJ  [_b  +  AloJ 

(4.30) 

and  the  independent  model  with  At(t),  A2”(t)  gives  for  0  <  rt  <  r2 
W  rl*r3  )s  p(  —  rl»  ^3  ^  r7  I  X7oi  ao(u)*  ai(u)*  b;  u  €  [0,oo)  ] 


o(Ti)+ai(ri) 


o(T3)+ai(Tj) 


(4.31) 
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It  is  easy  to  verify  that  theorem  4.4  holds. 

Theorem  4.4 

Let  F,nd(  rltr2  ),  Fcdf(  rlt  r2  )  ,  and  FLS(  rlf  r2  )  be  defined  as  above.  Then 

F inq(  ru  t2  )  <  Fcdf(  tx»  r3  )  <  F LS(  rlt  r2  ) 
for  rj  ,  r2  >  0  and  Ala,  A2e  >  0. 

Thus  the  bivariate  survival  function  of  2-component  parallel  red  undent  systems  increases  with 
the  degree  of  dependence  between  cumulative  EFFs  described  by  identically  distributed  gamma 
processes.  This  result  extends  to  m-component  systems  and  is  expected  based  on  the  results  of 
section  3.4  and  the  limiting  arguments  in  section  4.1.1. 
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APPENDIX  A 


a,  ar(t) 


<*.»  <*<(*)»  •=> . m 


aji  i= o.x.... 


a,  a(t) 


b 


0i  0 hi  k  — 0,1 . * 


CDF 


c*,  *=o,i . J 


NOTATION 


Shape  parameter  of  a  gamma  distribution  or 
process. 


Shape  parameter  of  a  gamma  distribution  or 
process,  referring  to  component  i. 

Shape  parameter  of  a  gamma  distribution  or 
process,  referring  to  the  EFF  rj;. 

Derivative  of  the  shape  parameter  «f  a  gamm*^ 
distribution  or  process. 

Scale  parameter  of  a  gamma  distribution  or 
process. 


Unknown  coefficient  or  known  distribution 
parameter. 

Refers  to  models  for  EFFs  ^(t)  and  jja(t)  that 
use  the  Cherian-David-Fix  bivariate  density 
(see  Chapter  2). 

Factor  used  in  section  3.1  to  describe  the  effect 
of  the  4th  stress  on  the  system. 
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o(t) 


Environmental  Factor  function  (EFF)  common  to 
all  components 


i)j,  j= EFF  over  time  interval  [  t^,  t^j  ),  common  to 

both  components  (see  Chapter  3). 

7,(t),  «= 1 . m  EFF  for  component  i 

G,  Gj,  Gt,  j— o.i,...,  A  distribution  function  for  a  random  variable. 

*=o,i . * 


G(  a,  b  ) 


Gpr(  <*(t),  0  ) 

Gpr(  a(t),  0{i)  ) 


h(t) 


H(t) 

«,  «  =  0,1, ...,m 


Gamma  distribution  with  density 


,  v  u®  v®— i  bx 
dG(x)  =  ^ — *T(a) - ’  X  >  °‘ 


Gamma  process  with  parameters  a(t)  and  0. 


Extended  gamma  process  with  parameters  a(t) 
and  0(t). 


Attenuation  function  for  the  shot-noise  model  of 
section  3.2. 

H(t)  =  |  h(u)  du. 

Subscript  referring  to  component  i. 


I*(t;),  *=o,i . *, 

i=o,i,... 


Indicator  variable  used  in  section  3.1. 

w  =  1  if  ith  stress  is  present  at  time  tj, 

0  otherwise. 


IND 


Refers  to  independent  models  for  EFFs  »7,(t)  and 
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h  j  =  0,1,... 


k,  k  =  0,1,... ,3 


Subscript  referring  to  time,  related  to  interval 
[  ti»  fci  +  i  )• 

Subscript  referring  to  a  stress  k  or  a  factor  ck 
used  to  describe  the  effect  of  stress  k. 


^«(^)i  -  — 1,2,. ..,m 

A(t) 

Aj(t),  «  =  1 . m 

\*  \*  \* 

'’ll  *12 


Baseline  failure  rate  for  component  L 

Total  failure  rate  for  component  i. 

A,(t)  =  Aio(t)  7i(t) 

Cumulative  failure  rate  for  a  single  component; 
A(t)  =  |  A0(u)  ?(u)  du. 

Cumulative  failure  rate  for  component  L 

Failure  rate  parameters  of  the  Marshall-OIlcin 
bivariate  exponential  distribution. 


LS 


Refers  to  model  derived  in  Lindley  h. 
Singpurwalla  [1986]. 


m,  m(t) 

M(t) 

nn  *4 

N j,  i= o,i,... 


Intensity  function  for  a  Poisson  process. 


M(t)  =  |  m(u)  du. 


Indices  related  to  times  rt  and  r},  such  that 
0  <  t»,  <  Tj  <  tnj+1  <  t«,  <  r7  <  tnj+l 


Number  ot  stresses  present  during  interval 
f  *>’  kj  +  i  )• 
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p»y*  q*>;  *=o.i . *» 

>  =  0,1.... 


is  the  probability  that  the  kth  stress  is 
present  during  time  interval  [  t^,  t  .^.j  ). 

«*>  =  1  -  p*>- 


_t  — ( tltt2,... ) 


Vector  of  (known)  transition  times  for  the 
piecewise  continuous  model  of  Chapter  3. 


I  li’  S  +  i  )' 


Time  interval  over  which  the  EFF  is  constant  in 
the  models  used  in  Chapter  3. 


T,-»  i 


Lifelength  of  component  i 


Time  at  which  a  univariate  distribution  is 
evaluated;  r  >  0. 


rU  r3 


Times  at  which  a  bivariate  or  marginal 
distribution  is  evaluated;  rlt  r7  >  0. 


At 

Y(t) 

Yj(t)  i  =  l . m 


At  =  [  tj,  tj  +  x  )  V  j>0  if  time  intervals  are  ^ 
uniform. 

Cumulative  EFF;  Y(t)  =  f  »7(u)  du. 

J  o 

Cumulative  EFF  for  component  ». 
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APPENDIX  B 


DETERMINISTIC  EFFs  WITH  UNKNOWN  COEFFICIENTS 


In  the  introduction,  we  stated  that  a  routine  extention  of  the  LS  model  is  to  describe  the 
dynamic  environment  using  an  environmental  factor  function  rj(r)  that  is  deterministic  with 
unknown  coefficients.  A  polynomial  function  of  time  is  used  to  demonstrate  general  results. 

Suppose  that  the  effect  of  a  dynamic  environment  on  the  system  over  a  finite  range 
[  0,  T’  )  can  be  captured  through  a  polynomial  environental  factor  function  rj{r).  We 
suppose  that  the  degree  of  the  polynomial,  p ,  is  known  or  assumed  but  the  coefficient  values 
are  unknown.  We  denote  the  coefficients  as  0k.  k  =  0,1,  .., p,  denote  the  vector  of 
coefficients  as  0  =  (  0Q,  0X,  ...,  0P  ),  and  describe  our  uncertainty  about  each  0k  with  a 
distribution  Gfc.  We  further  assume  that  the  0k's  are  mutually  independent  V  k, 
k  =  0,1 Thus  the  EFF  rj(r)  =  £  0k  r*.  P  >  0,  r  G  [  0,  T’  ). 

4=0 

Since  we  restricted  the  range  of  r  to  an  interval  [  0,  T’  ),  we  can  allow  negative 
coefficients;  thus  the  EFF  rj( t)  may  increase  or  decrease  within  the  interval  [  0,  T’  ).  The 
support  of  the  0ka  may  be  nonnegative  or  restricted  such  that  rj(r)  >  0  V  r  €  [  0,  T’  ). 
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B-l.  A  Polynomial  EFF  with  Gamma  Distributed  Coefficient* 


Suppose  that  our  uncertainty  about  each  of  our  coefficients  /?k  can  be  described  using  a 
gamma  distribution  with  known  parameters  ark  and  bk.  Then  our  system  survival  function 
will  be  for  any  r  €  [  0,  T’  ) 

P(  r  )  =  P[  T  >  t  |  Alal  A3o,  ak,  bk,  *=o,i . *  ]  = 


n 


*=0 


(t+l)bt 

r*  +  1  +  (*+l)bk 


«* 


+ 


p 

n 


(*+l)bt 

r‘  +  1  +  (i+l)bk 


<** 


p 

n 

*=0 


(*+l)bk 


(Al0+A3o)  r*  +  1  +  (fc+l)bk 


«* 


with  a  bivariate  survival  function  for  rlt  r3  6  [  0,  T’  ) 

P(  rii  rj  )  =  P[  Tj  >  rlt  Tj  >  r3  j  Ala,  A3o,  ak,  bk,  *=o,i . »  ] 


p 

n 

*=0 


'i«'  i 


rf  +  1  +  A2oi 


*  +  i 


+  (*+l)b» 


Our  marginal  survival  functions  will  be  for  ri  €  [  0,  T’  ) 
)  =  P[  T,  >  t{  |  Aio,  ok,  bk,  .  .  . . ] 


p 

n 

*=0 


(*+l)bk 


A, 


,r‘+l  +  (t+l)bk 


•  =  1,2 
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B-2.  A  Polynomial  EFF  with  Coefficients  Distributed 


as  a  Truncated  Normal 


Suppose  that  our  uncertainty  about  each  of  our  coefficients  0k  can  be  described  using  a 
truncated  normal  distribution  with  known  parameters  pk  and  <rk,  truncated  from  below  at 
some  value  0k.  We  denote  this  distribution  as  TN  „»(  nk,<r\  ).  Our  density  function  for 

Pk 

each  0k  will  be 


dGk(0k)  =  - ■—  exp  f" - ^  (0k  ~  ,  Pk  >  0k  , 

ck  L  J 


=  0, 


Pk<Pl> 


where 


oo 

=  [  -p— ;  «p  [--b  ( x  -  dx 
J  iinrl  L  J 

Pk 


We  need  to  choose  our  lower  truncation  point  Pk  of  the  distribution  Gt  in  such  a  way  as  to 
guarantee  that  »;(r)  >  0  V  r  6  [  0,  T’  ).  To  do  this,  we  establish  the  following  proposition: 


Proposition  B-l:  Define  t*  such  that  V  Pk(^*)k  =  m*n 

- — -  *=o  t€[0,T’] 

£  >  0,  then  rj(t)  >  0  V  t  >  0. 

*= o 

Proof  of  Proposition  B-l:  »j(t)  =  £  /M*  >  51  Pk^  ^  E  £*(**)*  >  0 

4=0  ~  k=0  4=0 


E  01 1*. 

4=0 


If 
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Thus,  to  insure  that  t)(t)  >  0,  we  can  adjust  the  values  0%  by  trial  and  error  until 

P 

12  Pk(t*)k  >  0-  Note  that  0*  must  be  positive  since  17(0)  >  0.  We  provide  several 

k—0 


examples  below  on  how  to  find  each  0\.  Obviously,  this  procedure  will  become  tedious  for 
large  values  of  p  unless  we  set  most,  if  not  all,  0\  >  0. 

Example  1:  p  =  1  (linear  function) 

We  know  that  0t  >  0.  If  0*  >  0,  then  £  0\rh  >0  V  r  >  0.  Suppose  0f  <  0.  By 

k—0 


definition,  0*  +  0*t*  =  min  0t  +  0ft.  A  linear  function  is  minimized  at  either 

t€[0,T’] 

endpoint;  if  0*  <  0,  then  the  function  is  decreasing  and  its  minimum  is  reached  at  T’.  Thus 
the  condition  is  met  for  any  0*  >  ^ 0  . 

Example  2:  p  =  2  (quadratic  function) 

Again,  we  know  that  0*  >  0.  If  0*  and  0\  >  0,  then  £  0l(t*)k  >  0  .  Suppose  that 

*=o 

0*  <  0  and  0$  <  0.  Then  0t  +  0*t  +  0*t2  is  minimized  at  t  =  T\  Therefore,  one  can 
select  any  0*  and  0*  such  that  0*  +  0%T'  >  •  Now  suppose  that  0*  <  0  and 

0*  >  0.  A  quadratic  will  be  minimized  at  the  endpoints  or  at  a  point  where  the  first 
derivative  equals  0.  Thus  t*  =  0,  T’,  or  -  0*/20l.  We  know  that  if  t*  =0,  then 
12  0*(4*)*  =  0*  >  0.  Any  value  of  0*  and  0\  such  that  £  £j(t*)‘  >  0  for  both 

*=0  k—0 

t*  =  T’  and  t*  =  —  0*/20\  will  be  satisfactory.  The  latter  condition  also  holds  if  we 


suppose  that  0*  >  0  and  0%  <  0. 


After  we  have  determined  the  parameters  0\  for  every  k,  k  =  0,1,. ...p,  then  we  can 
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easily  calculate  the  system  survival  function  for  any  r  >  0: 
F(  r  )  =  P[  T  >  r  |  A1)M  \3o,  ftk,  a2k,  0\,  *=o,i . i  ]  = 


n,0^)  +  A  g;<  ^ ,  +  fio  G>(  (Au+t+1>  -*+‘ ) 


f  A,  r*'H 

,n  M  Hfi- )  “p 


„  '*+1  ,  *  <T5T^’ 

^  k+ 1  +  - 


t+1 


n  Kk(  V?__ )  exp 

*=0 


-/i* 


l+,  -2  (*2®  r*+\2 

A,„r‘+1  .  **(-*+1  ) 


*2o 

T+T 


A  v  ,  (A10  +  A2o)  r‘+1 

fl0  Kt(  — i+i — 


•  exp 


..  (Aio  +  A2o)  r‘+l 

_/i‘  f+1 - 


2  /(*lo  +  *2o)  r*+\2 

( - E+I - ) 

2 


If  we  define  $  as  the  standard  normal  cdf,  then 


c* 


»(■) 


=  i  _ 


-(/jt-<r|s)1 

&h 


J’ 


and 


K*(>) 


1  -  * 

0k-(t*k-<r\») 

. _ Zl _ - 

1  -  * 

n-Hk 
L — — J 
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na* 


The  bivariate  survival  function  is,  for  0  <  rt  <  r2  <  T’: 

F(  ri»  r2  )  —  Tj  >  T|,  T2  >  r2  |  Alo,  A3o,  /ik,  <rt,  *=o,i . t  ]  = 


'i « 


*+i  ,  ,  r*+1 

rl  t*2<>  r3 


k+ 1 


) 


exp 


*+i 


~P* 


Mo 


ra 

1+1 


*+i  (rl  (*lf_ 
+  - 


*+i 


+ A 


3+r 


*+i 

2o  T7  \2 


with  marginal  survival  functions  for  r{  €  [  0,  T’  ) 

F,(  Ti  )  =  p[  T,-  >  r,  |  Aio,  nk,  a],  *=o,i . *  ]  = 


n  k*( 

*=o 


A-  r‘+1 

AtO  7  i 

4+1 


)  ««P 


“P* 


A io  T 


k+l 


*1  r' 


*+l 


TFT 


*+r 
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